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EXECUTIVE  SUMMARY 

Ice  is  extremely  variable,  and  because  waves  and  ice  may  interact  in  a  number  of  ways 
depending  on  the  ice  and  wave  characteristics,  wave  prediction  in  ice  covered  waters  is  a 
difficult,  and  in  some  ways  intractable,  problem.  A  review  of  wave-ice  interaction  is  made, 
and  special  attention  is  given  to  formulations  of  the  dispersion  relation.  The  dispersion 
relation  determines  the  wave  attenuation  for  non-conservative  dissipation  schemes  and 
contributes  to  a  change  of  wave  height  (and  direction)  analogous  to  shoaling  and 
refraction.  A  method  for  jointly  measuring  dispersion  and  attenuation  in  ice  is  proposed 
and  it  is  hoped  that  two  ongoing  ONR  Departmental  Research  Initiatives  will  provide  an 
opportunity  to  test  and  evaluate  models  of  dispersion  in  ice.  In  terms  wave  effects  on  ice, 
a  number  of  new  primitive  ice  models  have  been  recently  developed,  these  are 
summarized,  but  much  research  is  needed  towards  understanding  the  mechanical 
properties  of  ice.  Also,  additional  analysis  of  wave  attenuation  during  a  high  wave/ice 
breakup  is  presented  as  a  companion  for  Collins  et  al.  [2015], 
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1  INTRODUCTION 

This  document  is  concerned  with  the  interaction  of  wind-generated  sea  surface  gravity 
waves  (henceforth  waves)  with  sea  ice.  The  physics  of  wave-ice  interaction  is  equally 
applicable  in  both  of  the  Polar  Regions,  although  the  oceanography  of  the  Arctic,  the 
exclusive  location  of  the  recent  Office  of  Naval  Research  (ONR)  sponsored  field  work  and 
the  focus  of  this  report,  is  quite  distinct  from  the  Antarctic.  In  simple  terms,  the  Arctic  is 
ocean  surrounded  by  land  and  the  Antarctic  is  land  surrounded  by  ocean.  This  difference 
in  geography  plays  a  role  in  ice  extent  as  a  result  of  wave  action:  the  Arctic  is  mostly  ice- 
covered  in  the  winter  so  fetch,  and  hence  wave  action,  develops  primarily  in  warmer 
summer  months.  Waves  may  facture  the  ice  which  increases  the  lateral  surface  area, 
increased  surface  area  accelerates  melt  for  smaller  floes.  An  increased  fraction  of  open 
water  also  leads  to  albedo  feedback  effects  and  increased  fetch.  The  increased  fetch 
further  increases  the  potential  for  the  development  of  more  damaging  waves.  Overall, 
waves  work  to  reduce  the  ice  in  the  Arctic,  while  in  the  Antarctic  the  most  energetic 
waves  occur  in  the  winter  at  which  point  ice  fracturing  may  contribute  to  increasing  the 
ice  extent.  This  is  because  waves  are  facilitating  the  outward  export  of  ice  as  broken  floes, 
with  the  leads  between  floes  subsequently  freezing  over  during  these  cold  months  [Squire 
and  Montiel,  2015], 

ONR  has  sponsored  two  related  Department  Research  Initiatives  (DRI):  the  Marginal  Ice 
Zone  Program  and  the  Sea  State  and  Boundary  Layer  Physics  of  the  Emerging  Arctic  Ocean 
(Sea  State)  [Thomson  et  al.,  2013],  The  purpose  of  this  document  is  to  review  the  current 
body  of  science  in  anticipation  of  new  results  from  these  DRIs  as  well  as  provide 
documentation  of  additional  analysis  of  a  dataset  from  the  University  Center  in  Svalbard. 
This  document  is  necessarily  limited  in  scope,  for  more  comprehensive  reviews  of  wave- 
ice  interaction  in  the  MIZ  please  see  chapter  6  of  Wadhams  [2000],  Squire  et  al.  [1995], 
and  Squire  [2007],  Parts  of  section  2.2  have  been  expanded  in  an  independent  manuscript 
in  progress  of  publication  [Collins  et  al.,  2016], 

1.1  THE  MARGINAL  ICE  ZONE 

Large  regions  of  the  Arctic,  particularly  in  the  winter  months,  experience  total,  compact 
ice  cover  (so-called  pack  ice).  In  between  the  ice-free,  open  seas  and  the  pack  ice  exists  a 
gradient  of  ice.  This  transition  zone  is  the  marginal  ice  zone  (MIZ).  The  MIZ  is 
characterized  by  ice  which  weakens  in  concentration  towards  the  open  sea,  often  with 
decreasing  floe  sizes  and  thicknesses.  As  the  ice  concentration  increases,  so  does  the  ice 
effect  on  wave  propagation.  Deep  in  the  ice  pack,  wave  energy  is  severely  attenuated  and 
therefore  wave  effects  on  ice  are  limited. 

As  is  easily  observable  with  satellite  remote  sensors  (which  have  been  operational  since 
1978),  on  the  scales  of  seasons,  the  MIZ  retreats  into  the  interior  of  the  Arctic  as  ice  melts 
during  the  summer  months  reaching  a  minimum  on  average  in  September.  The  ice  extent 


Manuscript  approved  February  10,  2016. 


1 


returns  during  the  winter  as  the  Arctic  water  refreezes.  The  ice  that  survived  from  the 
previous  year  is  known  as  "multi-year  ice"  while  new  ice  is  called  "first-year  ice".  This 
distinction  is  important  because  the  mechanical  and  thermodynamic  properties  are  quite 
different  between  the  two.  On  the  scale  of  climate,  observations  show  that  more  and 
more  multi-year  ice  is  melting  (implying  decreasing  minimum  ice  extent),  and  thus  the 
ratio  of  first-year  ice  to  multi-year  ice  is  on  the  rise.  This  has  led  to  the  speculation  of  a 
possible  ice-free  Arctic  in  the  future,  figuratively  and  literally  opening  up  a  new  ocean  to 
the  world  with  "new"  natural  resources  available  for  exploitation.  Already,  there  is  more 
navigable  Arctic  Ocean  due  to  the  historically  low  summer  ice  extents. 

Therefore,  there  is  an  impetus  to  better  understand  the  retreat  of  sea  ice  including 
possible  connections  of  scales  through  feedback  loops.  Besides  the  large  scales  previously 
mentioned,  wave-ice  interaction  occurs  on  mesoscales  (e.g.  MIZ  migration  during  large 
storms),  sub-mesoscales  (smaller  storms  and  gales,  changes  in  floe  size  distribution),  and 
micro-scales  (energy  fluxes,  fracture  events).  Little  is  known  about  the  smaller,  sub-meso- 
and  micro-scales  because  of  lack  of  comprehensive  observations.  The  main  idea  behind 
the  Arctic  Sea  State  DRI  is  to  observe  at  these  scales  and  integrate  the  findings  into 
current  operational  oceanography. 

1.1.1  Wind-Generated  surface  gravity  waves 

Waves  are  a  ubiquitous  and  fundamental  part  of  nature.  All  bodies  of  water  on  Earth 
mediate  the  propagation  of  surface  waves,  and  understanding  of  their  physics,  dynamics, 
and  interactions  is  an  essential  prerequisite  for  describing  the  Earth  environment  and 
climate  system.  For  many  applications,  particularly  sea-keeping  and  ocean  operations, 
medium-range  forecasts  with  hourly  resolution  are  sufficient.  At  these  scales,  information 
about  the  waves  is  reported  in  terms  of  the  wave  spectrum,  i.e.  individual  waves  are 
treated  statistically  not  deterministically,  though  the  spectrum  itself  evolves 
deterministically.  So-called  third  generation  (3G)  spectral  wave  models  were  developed 
for  this  forecasting  of  this  class,  and  their  operational  implementation  have  become  one 
of  the  great  success  stories  of  modern  geoscience. 

1.1. 1.1  3G  WAVE  MODELS 

All  3G  wave  models  (e.g.  SWAN  "Simulating  WAves  Nearshore"  [Booij  et  al.,  1999],  WW3 
"WAVEWATCH  III"  [Tolman,  1991;  Tolman  and  the  WAVEWATCH  III®  Development  Group, 
2014])  are  based  on  the  phase-averaged,  action  balance  (also  known  as  radiative  transfer) 
equation: 

DiV  v  5f 

Dt  L-i  a 
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This  describes  the  conservation  of  wave  action  (spectral  density  divided  by  frequency) 
where  the  material  derivate,  ^  -  change  in  time  and  divergence,  balances  the  sum  total  of 
sources  and  sinks  of  wave  energy.  Si,  and  a  is  the  intrinsic,  radial  wave  frequency. 

On  the  open  ocean  and  in  deep  water,  the  main  source  terms  are  input  of  energy  by  wind 
(Sin),  dissipation  of  energy  by  breaking/white  capping  (Sds),  and  redistribution  of  wave 
energy  through  weakly  non-linear,  wave-wave  interactions  (Sni)  (e.g.  Cavaleri  et  al.,  2007). 
Recent  work  has  focused  on  refinement  of  Sin  and  Sds  (see  Ardhuin  et  al.,  2009,  Rogers  et 
al.,  2012,  and  references  within),  and  more  realistic  implementation  of  Sni  (e.g.  Tolman 
2013;  Rogers  and  Van  Vledder,  2013). 

When  waves  propagate  into  the  MIZ,  there  is  strong  a  modification  of  these  source  terms 
(though  the  details  are  not  known).  Besides  modification  of  open  water  source  terms,  the 
wave  interaction  with  sea  ice  necessitates  a  new  source  term,  SiCe.  The  first 
implementation  of  ice  effects  on  waves  was  a  simple  partial  blocking  of  energy  flux  scaled 
by  ice  concentration  [Tolman,  2003],  rather  than  via  physics  (S\ce).  This  was  recently 
improved  upon  by  the  implementation  of  physics-based  S\ce  schemes  [Doble  and  Bidlot, 
2013;  Rogers  and  Orzech,  2013;  Rogers  and  Zieger,  2014], 

1.1.2  Sea  Ice 

Though  somewhat  variable,  the  salinity  of  oceanic  waters  differentiates  it  from  its  fresh 
water  counterpart  by  35  "salt"  grams  per  liter.  These  salts  fundamentally  alter  the 
thermodynamics  of  the  water  including  the  freezing/melting  point.  At  this  typical  salinity 
(and  at  the  pressure  of  1  atmosphere),  the  freezing  point  of  sea  water  is  -2°  C. 

The  mechanical  differences  between  first-  and  multi-  year  ice  were  previously  mentioned. 
When  sea  ice  forms,  the  water  freezes  and  the  salts  are  concentrated  into  brine  channels. 
These  brine  channels  add  to  the  porosity  of  ice.  It  turns  out  that  porosity  is  inversely 
proportional  to  the  flexural  strength,  i.e.  very  porous  ice  has  weak  flexural  strength.  As  ice 
warms  in  the  spring,  the  network  of  brine  channels  may  connect  and  drain  out  of  the  ice. 

If  the  ice  survives  the  melt  season,  then  it  is  much  less  porous  after  refreezing.  In  this  way 
multiyear  ice  possess  much  more  flexural  strength  and  therefore  is  more  resistant  to 
breaking  [Timco  and  Weeks,  2010],  Much  of  the  difficulty  in  determining  the  appropriate 
wave-ice  interaction  is  a  result  in  uncertainty  in  mechanical  ice  properties.  These 
properties  are  difficult  (and  in  some  cases  impossible)  to  measure  in  the  field.  Various 
mechanical  properties  are  important  to  the  study  of  ice  -  wave  interaction  including  (but 
not  limited  to):  thickness,  density,  porosity,  salinity,  temperature.  Young's  modulus,  shear 
modulus,  Poisson's  ratio,  creep  and  fatigue  properties,  and  breaking  strength.  Due  to 
difficulties  of  in  situ  measurements  of  these  properties,  these  are  better  known  for  first 
year  ice  than  for  multiyear  ice.  Even  in  a  small  area  of  first  year  ice,  the  mechanical 
properties  may  be  subject  to  significant  additional  dependencies  in  vertical  axis  due  to 
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change  in  temperature  and  salinity  or  in  time  due  to  creep,  fatigue,  and  eventual  failure 
[Squire  and  Montiel,  2015], 


In  addition  to  first-  and  multi-year  ice  solid  pack  ice,  there  are  many  more  forms  of  sea  ice 
forged  by  interactions  with  the  environment  and  all  of  these  may  be  varied  and  nuanced 
in  nature.  The  fact  that  the  Inipaiq  Eskimos  of  Wales,  Alaska  have  120  words  for  sea  ice 
[Krupnik  and  Weyapuk  Jr,  2010]  reflects  the  many  variations  of  its  form  in  nature.  For  a 
review  of  the  properties  of  sea  ice  in  the  Arctic,  see  Wadhams  [1981],  and  an  engineering 
context  see  Timco  and  Weeks  [2010], 

Since  ice  is  so  varied,  it  is  useful  to  define  homogeneity  of  ice  conditions  in  two  ways.  If, 
on  the  large  scale,  a  particular  mechanism  of  wave-ice  interaction  is  dominant,  but  the  ice 
characteristics  are  varied,  then  this  may  be  thought  of  as  a  qualitative  large-scale 
homogeneity.  An  example  would  be  a  large,  relative  thin  ice  sheet  with  a  thickness  which 
slowly  changes  over  space.  Small-scale  homogeneity,  by  contrast,  is  more  strictly  defined. 
In  the  previous  example,  strict  (small-scale)  homogeneity  would  imply  uniform  ice 
thickness.  If  multiple  regimes  of  wave-ice  interaction  are  appropriate  within  the  scale  of 
interest  (e.g.  multiple  scattering  by  individual  floes  abruptly  transitioning  to  solid  pack 
ice),  this  violates  any  sense  of  homogeneity.  Obviously,  the  appropriateness  of  the 
homogeneity  assumption  will  depend  on  the  scales  of  interest. 

1.2  HISTORICAL  NOTE 

Man  has  explored  frozen  waters  near  or  in  the  Arctic  since  the  time  of  the  Viking  Era.  In  a 
quest  for  shorter  shipping  routes,  there  were  several  expeditions  into  the  Arctic  Ocean  in 
the  1800s.  Ships  would  become  stuck  in  the  sea  ice,  some  relenting  under  the  immense, 
crushing  pressure.  The  Fram  was  specially  built  for  withstanding  this  pressure;  to  survive 
the  grip  of  the  ice  and  drift  along  within  it.  During  its  maiden  voyage  (1893-1896)  the  crew 
sailed  deep  into  the  ice  with  the  intentions  of  putting  Fram  to  this  test.  During  this  time 
the  crew  managed  many  to  record  many  scientific  soundings  (e.g.  pressure,  depth, 
location).  The  first  explorers  to  reach  the  North  and  South  Pole  achieved  the  feat  in  1908 
(or  possibly  1909)  and  1911,  respectively. 

By  the  mid-1940s,  the  Soviets  were  by  far  the  leaders  in  Arctic  science,  having  had  an 
almost  continuous  scientific  presence  since  1935.  Professor  N.  N.  Zubov's  book  on  Arctic 
ice  (which  was  eventually  translated  into  English)  was  perhaps  the  first  treatise  on  the 
subject  [Zubov,  1945], 

Many  passages  in  Zubov's  book  are  remarkable;  here  we  reproduce  several  which  are 
directly  related  to  the  interaction  of  waves  and  ice: 

The  deflection  of  ice  under  the  weight  of  a  load  is  ordinarily  compared  to 
the  problem  of  deflection  of  an  elastic  plate  on  an  elastic  foundation.  But 
this  problem  presupposes  uniformity  of  the  plate  and  ice  is  extremely 
heterogeneous,  both  vertically  and  horizontally...  Finally,  when  loads  are 
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moved  on  ice,  the  ice  deflects  in  conformity  with  the  "water"  wave  which 
forms  under  it.  All  these  facts  taken  together  create  unusual  difficulties 
for  a  mathematical  analysis  and  make  it  necessary  to  have  recourse  to 
formulas  which  even  though  approximate,  still  satisfy  the  practical 
requirements. 

His  description  of  the  state  of  affairs  for  the  treatment  of  the  waves  propagating  through 
ice  is  astonishingly  appropriate  for  modern  wave-ice  interaction  modeling.  Since  the  time 
of  Zubov's  book,  many  analytical  problems  have  been  solved  with  sophisticated 
treatments  of  waves  propagating  through  non-uniform  ice  (e.g.  Squire,  2007).  However 
the  complex  mathematical  treatments  are  not  appropriate  for  application  in  phase- 
averaged,  spectral  wave  models  which  in  practice  are  the  workhorses  responsible  for 
forecasting  waves  in  the  Arctic  and  thus  many  simplifications  are  necessary. 

An  anecdote,  related  by  Zubov  via  an  earlier  publication  by  Bernstein,  might  be  the  first 
recorded  observations  of  waves  under  ice  cover: 

Observer  repeatedly  noted  the  appearance  of  wind  ripples  on  the  thin 
Ice.  Bernstein  points  out  that  the  instrument  observations  carried  out  in 
1927  on  the  Volga  established,  without  a  doubt,  the  origin  of  the  wind 
fluctuations  of  the  ice  level.  Thus,  with  a  wind  of  13  m/sec.  (which 
corresponds  to  a  pressure  of  21  kg/m2  on  a  surface  perpendicular  to  the 
wind),  a  constant  agitation  of  the  Ice  (Irregular  periods  of  20  to  180  sec.) 
was  observed.  The  amplitude  reached  3  mm  (figure  86).  In  calm  weather, 
fluctuations  of  an  equal  order  were  not  observed. 

The  Volga,  running  through  central  Russia,  is  the  longest  and  largest  (by  discharge)  river  in 
Europe.  Observations  of  Arctic  ice  vibration  thought  to  be  connected  to  wind  disturbances 
were  published  by  Fakidov  [1934]  where  it  was  noted  that  vibrations  would  often  precede 
the  oncoming  wind.  This  might  be  explained  by  the  propagation  of  flexural-gravity  waves 
which  may  outrun  the  storm  which  generated  them. 

The  Soviet  Icebreaker  Gerogy  Sedov  was  sent  to  the  Arctic  in  1929  for  a  high  latitude 
expedition.  In  1937,  the  Sedov  became  stuck  in  the  ice  and  was  eventually  converted  into 
a  research  station  (the  first  North  Pole  research  station)  which  spent  812  days  drifting  in 
the  ice.  The  captain's  log  was  described  in  Zubov's  book  with  several  fascinating  entries 
which  may  be  the  first  accounts  of  wave  causing  the  breakup  of  Arctic  ice: 

On  26  January  [1938],  a  storm  set  in  which  lasted  six  days.  The  ice  field 
began  to  undergo  more  powerful  oscillations.  The  period  of  these 
oscillations  was  10  to  12  seconds  and  was  the  sum  of  the  waves  [sic] 
period  and  the  period  of  its  own  fluctuations.  The  inclination  of  the  ice 
field  reached  60  angular  seconds  or  more.  As  a  result  of  these 
oscillations,  tensions  were  generated  in  the  ice  field,  which  finally  caused 
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it  to  break  up  on  1  February  along  lines  approximately  perpendicular  to 
the  wind  direction.  A  large  surge  undoubtedly  caused  these  fluctuations 
and  thus  the  break-up  of  the  ice  field.  The  surge  was  caused  by  stormy 
winds  in  the  neighboring  ice-free  areas  of  the  Greenland  Sea  and  spread 
in  all  directions  according  to  a  general  law. 

Again,  according  to  Zubov: 

On  12  January  1940,  it  was  noted  in  the  log  of  the  Sedov: 

'2100,  significant  surge  30  m  -  when  the  ridge  passed,  ice  was  raised  a 
little  and  broke  open.  When  the  base  of  the  wave  had  passed,  the  ice 
rejoined  with  a  creak,  characteristic  of  jamming.  The  phenomenon 
occurred  with  a  periodicity,  characteristic  for  a  heavy  sea  of  9-10  sec.' 

Over  the  years,  seismic  or  "gravity  meter"  observations  followed  which  correlated 
millimeter  (or  less)  movement  of  the  ice  with  wind  or  storms  [Crary  et  al.,  1952;  Hunkins, 
1962;  LeSchack  and  Haubrich,  1964], 

If  the  ice  floes  are  much  smaller  than  the  wave  length  of  the  impending  ocean  wave,  then 
gravity  is  still  the  sole  restoring  force  (though  there  is  additional  inertia  due  to  the  mass  of 
ice)  and  these  waves  are  essentially  ocean  waves  modified  by  the  presence  of  ice.  As 
water  waves  make  their  wave  into  much  larger  floes  (than  the  wavelength)  or  solid  pack 
ice,  where  the  model  of  a  large  elastic  sheet  or  beam  may  be  appropriate,  then  they 
prorogate  through  the  ice  as  coupled  wave  modes  known  as  flexural-gravity  waves.  These 
flexural-gravity  waves  may  be  thought  of  as  "ice  waves"  as  now  the  restoring  force  is 
gravity  plus  the  flexural  properties  of  the  ice  sheet  or  beam,  in  other  words,  ice  is  the 
interfacial  medium  which  supports  the  waves,  and  it  was  these  ice  waves  which  were  first 
measured.  The  first  account  of  ocean  waves  propagating  through  the  marginal  ice  zone 
and  ice  fields  comes  from  Robin  [1963],  Studies  postdating  these  will  be  referenced  in  the 
review  below  (section  2)  and  cataloged  in  the  context  of  wave  height  in  Appendix  B. 
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2  REVIEW 


2.1  ICE  EFFECTS  ON  WAVES 

As  ice  is  varied  in  nature,  the  effect  of  ice  on  waves  is  also  varied  (e.g.  Campbell  et  al., 
[2014]).  The  most  documented  property  is  attenuation.  Attenuation  is  the  result  of  many 
individual  mechanisms  responsible  for  dissipating  wave  energy  as  waves  interact  with  ice. 
These  mechanisms  are  typically  split  into  two  broad  categories:  1)  conservative  (SiCe,c)  and 
(2)  non-conservative  attenuation  mechanisms  (S\ce,nc)  [Rogers  and  Zieger,  2014],  In  the 
context  of  a  spectral  wave  model,  SiCe,c  includes  wave  scattering  and  reflection,  so 
although  wave  energy  is  attenuated  along  the  main  propagation  direction  into  the  ice,  the 
energy  is  being  redirected  (e.g.  backwards  toward  the  open  water)  rather  than  dissipated. 
For  Sice.nc,  waves  also  attenuate  along  the  main  propagation  direction,  but  the  wave 
energy  is  essentially  lost  to  unresolved  dissipation  processes.  Some  losses  are  caused  by 
differences  between  local  horizontal  motion  of  the  ice  and  water,  implying  a  frictional 
effect  which  goes  to  turbulence  and  then  heat  in  the  water.  In  the  case  of  ice  floes,  this 
friction  results  in  horizontal  motion  of  the  ice  floes,  and  if  these  floes  are  colliding  or 
ridging,  wave  energy  will  be  lost  to  friction  and  even  noise  and  potential  energy.  Normal 
pressures  at  the  water/ice  interface,  on  the  other  hand  do  not  inherently  dissipate  wave 
energy,  but  if  there  are  energy  losses  in  the  ice  caused  by  the  vertical  motion  of  the  ice 
(e.g.  hysteresis  or  internal  friction),  then  energy  is  again  lost  to  heat1. 

In  general  there  are  two  overarching  regimes  of  interaction.  One  is  when  the  horizontal 
extent  of  the  ice  floe  is  large  compared  to  the  typical  wavelength.  In  this  limit,  one  may 
model  ice  as  an  elastic  sheet,  then  wave  may  propagate  through  the  ice  as  flexural-gravity 
waves.  Dissipation  in  this  regime  may  be  due  to  formation  of  turbulent  boundary  layers 
beneath  the  ice.  The  second  is  when  the  horizontal  extent  of  the  ice  floe  is  small 
compared  to  the  typical  wavelength.  Then  the  waves  are  essentially  water  waves  with 
extra  mass  loading.  Grease  ice  may  damp  small  waves  via  viscous  interaction  of  ice 
particles.  Snell's  law  may  be  applied  to  understand  change  in  propagation  direction  and 
wave  height  of  the  incoming  waves  in  analogy  to  refraction,  and  similar  for  group  velocity 
in  analogy  to  shoaling  (more  on  this  in  section  2.1.4). 

2.1.1  Interaction  Mechanisms 

2. 1.1.1  CONSERVATIVE  ATTENUATION:  SCATTERING 

Scattering  is  a  universal  wave  process  by  which  waves  are  deflected  by  some  object.  A 
classic  example  is  the  Rutherford's  gold  foil  experiment,  where  he  discovered  the  atomic 
nucleus  from  scattering  patterns  of  alpha  particles  aimed  at  the  foil.  Similarly  ocean  waves 
may  be  scattered  by  obstacles,  including  ice.  Scattering  sends  part  of  the  wave  energy 


1  One  might  speculate  that  this  production  of  heat  in  the  ice  could  contribute  in  some  small  way  to  melting. 
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away  from  the  main  direction  of  propagation,  so  that  along  the  main  direction,  wave 
energy  is  attenuated.  If  ice  is  much  larger  in  scale  then  the  incident  wavelength,  then 
scattering  will  important. 

A  simple  ice  model  is  a  semi-infinite  thin  elastic  sheet  (so  called  Euler-Bernoulli  thin  elastic 
plate  [Squire,  2007]).  Waves  will  scatter  at  the  interface,  so  that  some  of  the  energy  is 
transmitted  into  the  ice  and  some  of  the  energy  is  reflected  back  out  to  sea.  As  waves 
propagate  through  realistic  ice,  there  is  scattering  due  to  imperfections  (i.e.  ridges  and 
sails)  and  discontinuities  (i.e.  change  in  ice  thickness,  cracks,  etc.).  Scattering  also  occurs 
upon  entering  and  leaving  leads  and  polynya  (open  water  regions  in  pack  ice  cover). 
Scattering  may  also  occur  in  the  interaction  with  a  single  floe,  and  many  of  these  single 
floe  interactions  can  be  synthesized  together  in  a  sophisticated  way  and  treated  as 
multiple  scatterers.  Many  complicated  mathematical  models  exist  for  scattering  and  the 
specific  problems  mentioned  above,  see  [Squire,  2007]  for  a  review. 

Typically,  to  solve  the  scattering  problem,  one  needs  detailed  information  about  the  ice 
floes  and  wave  phases  which  is  never  available  in  practice,  so  most  of  this  work  has  not 
been  translated  into  a  form  for  practical  application  in  3G  wave  models.  The  recent 
development  of  primitive  wave-ice  interaction  models  [section  2.2.4]  has  allowed  the 
scattering  attenuation  based  on  the  number  of  floes  in  the  domain  [Kohout  and  Meylan, 
2008], 

2. 1.1.2  NON-CONSERVATIVE  ATTENUATION:  DISSIPATION 

Alternative  to  or  in  concert  with  scattering,  one  may  consider  non-conservative 
dissipation.  The  main  mechanisms  for  non-conservative  dissipation  are  covered  in  Zhao  et 
al.  [2015]: 

1.  One  way  to  introduce  dissipation  is  through  viscosity  or  internal  friction 
(hysteresis)  of  the  ice.  For  example  modeling  the  ice  layer  as  a  continuum  with 
some  finite  viscosity,  that  is,  there  is  a  finite  ice  layer  has  some  rheological 
properties  would  be  appropriate  for  frazil  type  ice.  A  purely  elastic  ice  layer  does 
not  permit  this  type  of  dissipation.  However  viscosity  is  introduced,  it  is  a  two- 
layer  fluid  problem,  and  the  viscosity  of  the  ice  layer  or  water  layer  will  leach  wave 
energy.  Several  rheological  models  are  available  in  the  literature  (e.g.  Robinson 
and  Palmer  [1990],  Fox  and  Squire  [1994],  and  Keller  [1998]) 

2.  Another  way  to  dissipate  energy  is  in  the  viscous  sublayer  below  ice.  This  type  of 
dissipation  is  possible  even  when  assuming  that  the  ice  is  purely  elastic  (i.e.  no 
viscosity  in  the  ice).  Imagine  a  very  large  ice  sheet,  such  that  the  ice  does  not  move 
horizontally  to  match  the  wave  induced  horizontal  velocity.  The  interaction  of  the 
horizontal  water  motion  with  a  stationary  boundary  induces  a  turbulent  boundary 
layer  flow  directly  beneath  the  ice.  The  velocity  of  the  particles  goes  must  to  zero 
at  some  point,  depending  on  the  roughness  length  associated  with  the  underside 
of  the  ice.  Energy  is  dissipated  in  the  viscous  sublayer  close  to  the  ice.  This 


"turbulence  under  ice"  concept  has  been  used  in  the  sink  term  formulated  by  Liu 
and  Mollo-Christensen  [1988]  and  Liu  et  al.  [1991],  It  should  be  kept  in  mind  that 
their  derivation  is  primarily  to  predict  the  dispersion  relation  for  a  case  of  mass¬ 
loading  with  an  elastic  sheet,  and  the  "turbulence  under  ice"  scheme  was  added  as 
an  afterthought  via  an  eddy  viscosity  term  in  order  to  produce  dissipation.  It  is  also 
a  concept  used  in  the  physical  argument  of  De  Carolis  and  Desiderio  [2002],  who 
again  apply  it  as  if  the  turbulence  is  a  viscosity  in  the  water  layer. 

3.  Consider  again  an  ice  floe  which  is  large  compared  to  the  wavelengths.  As  the 
wave  propagates  as  a  flexural-gravity  wave,  the  ice  under  goes  compression  and 
decompression  twice  per  wavelength.  This  ice  will  react  with  some  inelasticity, 
meaning  the  bending  motion  alters  the  mechanical  structure  of  the  ice.  This  is 
known  as  elastic  hysteresis.  Related  to  this  a  mechanism  for  structural  failure  and 
eventual  breakup  of  ice.  At  some  number  of  cyclical  loads,  N,  the  ice  may  facture 
at  stresses  lower  than  the  yield  stress,  a  phenomenon  known  as  fatigue.  Another 
related  problem  is  the  permanent  strain  due  to  sustained  stress,  known  as  creep. 
Creep  is  an  issue  closely  related  to  temperature.  Wadhams  [1973]  provides  a 
formulation  of  the  creep  problem.  Another  way  to  think  of  these  processes  is  in 
terms  of  failure  mechanics.  Creep  is  ductile  or  plastic  failure,  where  the 
deformation  or  stretching  of  the  material  leads  to  failure  over  time.  Fatigue  is  the 
propagation  of  small  cracks,  where  the  crack  length  lengthens  over  time  and 
eventual  lead  to  failure.  The  creep  physics  used  by  Wadhams  [1973]  were  taken 
from  existing  knowledge  of  slow  movement  of  glaciers,  so  its  validity  for  the  rapid 
cyclical  motion  of  the  wave-in-ice  problem  is  not  known.  On  the  other  hand,  the 
ability  of  this  motion  to  cause  fatigue  in  the  ice  is  intuitive. 

4.  Now  consider  ice  floes  which  are  small  compared  to  the  wavelengths  such  that  the 
ice  floes  more  or  less  respond  to  the  horizontal  particle  motion.  If  the  floes  are 
densely  packed  then  there  will  be  collisions  between  floes.  There  are  a  number  of 
processes  resulting  from  inelastic  collision:  the  edges  of  floes  are  floes  may  be 
smashed,  grated,  pulverized  resulting  in  a  slurry,  or  the  floe  may  crack,  break,  or 
raft  (pile  on  top  of  each  other)  etc.  The  permanent  changes  in  ice,  including 
potential  energy  via  rafting,  must  drain  some  energy  from  the  wave  system.  Also, 
significant  noise  is  may  be  generated  as  floes  collide,  effectively  transporting  away 
some  energy  from  the  wave  system.  These  effects  are  of  such  complexity  that  they 
may  remain  strictly  in  the  realm  of  phase-resolved  and  floe-resolved  modeling,  to 
be  included  in  WW3-type  models  only  in  gross  and  highly  parameterized 
representations. 

One  may  take  lump  all  of  these  effects  together,  in  analogy  to  turbulence  closure,  into  an 
effective  phenomenological  viscosity  [Liu  and  Mollo-Christensen,  1988;  Wang  and  Shen, 
2010],  The  latter  intended  their  model  as  a  continuous  (i.e.  without  switching  between 
models)  representation  of  the  ice  cover  via  two  variables  which  are  specified 
independently: 
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The  viscosity  property  comes  from  the  frazil  ice  or  ice  floes  much  smaller 
than  wavelength.  Interaction  of  these  small  "particles"  and  their 
hydrodynamic  interaction  with  the  surrounding  water  create  an  effective 
viscosity  for  the  ice  layer.  The  elasticity  property  comes  from  the  rigidity 
of  ice  floes  in  which  floe  sizes  are  relatively  large  compared  to  the 
wavelength.  When  the  ice  is  consisted  of  frazil  or  small  ice  floes,  a 
viscous  parameterization  should  be  appropriate  for  the  ice  cover.  When 
the  ice  cover  is  a  continuous  ice  sheet,  an  elastic  parameterization  is 
appropriate. 


2.1.2  Attenuation 

Introducing  dissipative  processes  into  the  formulation  of  the  interaction  problem  may 
yield  direct  estimates  of  attenuation.  The  difficulty  is  that  various  dissipation  mechanisms 
(some  with  well-defined  models,  others  without)  may  or  may  not  be  appropriate  for  given 
conditions.  Some  require  homogenous  (even  isotropic)  conditions,  where  the  MIZ  is 
heterogeneous  in  general.  Many  mechanisms  may  be  active  at  the  same  time.  Apart  from 
dissipation,  wind  input  is  still  active  over  ice,  although  it  is  thought  to  be  drastically 
reduced  [Bidlot  et  al.,  2014],  There  may  be  normal  wave  generation  in  the  small  fetches  of 
leads  and  polynya.  The  sum  total  of  all  of  these  effects  are  rolled  up  into  the  calculation  of 
bulk  attenuation  of  wave  energy  typically  as  a  function  of  distance  into  the  ice. 


Figure  1  Image  of  the  MIZ  with  hypothetical  measurement  points  at  various  lengths  starting  at  the  open 
sea  edge. 


The  solutions  to  the  scattering  problem  as  well  as  measurements  indicate  wave  energy 
attenuations  exponentially  of  as  a  function  of  distance  [Wadhams  et  al.,  1988]: 


dSx 

dx 


=  ~aS  o 


(1) 


Here  So  is  the  wave  frequency  spectrum  at  position  x  =  0  and  a  is  the  attenuation 
coefficient.  This  is  solved  as  follows: 
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sx(f)  =  s0(f )  exp(-ax) 


(2) 


Imagine  the  situation  as  pictured  in  Figure  1.  There  are  three  measurement  locations. 
Starting  at  the  open  sea  edge  at  the  far  left  you  measure  in  the  incoming  wave  spectrum, 
So(f).  You  also  have  measurements  of  the  wave  spectrum  at  two  locations  deeper  into  the 
MIZ,  Sx(/).  In  this  case,  Eq.  2  may  be  inverted  to  yield  the  coefficient  of  attenuation: 

log(S* /S0) 


Scattering  models  and  non-conservative  dissipation  models  may  produce  a  direct 
estimate  of  the  attenuation  which  may  be  compared  to  these  measurements.  However, 
measuring  the  attenuation  in  this  way  does  essentially  provides  no  information  on  details 
of  dissipation  mechanisms  (hence  the  term  "bulk  property"  is  sometimes  used)  and  other 
sources  and  sinks  or  wave  energy  unrelated  to,  but  perhaps  strongly  modified  by,  ice  may 
be  active.  So,  it  may  come  to  no  surprise  that  it  is  not  a  trivial  task  to  compare  model 
estimates  of  a  with  measurements  nor  should  measurements  be  expected  to  agree  under 
different  conditions. 

The  ideal  measurement  scenario  is  to  have  an  unforced  swell  system  that  is  propagating 
directly  into  the  ice,  normal  to  a  clearly  defined  ice  edge  (or  floe  field  boundary),  with  in 
situ  measurements  at  several  points  inside  the  ice,  and  one  point  outside.  The 
measurements  should  be  directional  in  order  to  distinguish  the  effects  of  scattering  vs. 
dissipation.  Of  course,  a  real,  in  situ  measurement  effort  will  be  faced  with  wave  and  ice 
conditions  that  are  considerably  less  "tidy",  and  so  there  usually  must  be  some  educated 
guesswork  and  approximation. 

2.1.3  Dispersion  Relationship  in  Ice 

Attenuation  due  to  non-conservative  mechanisms  is  inseparable  from  dispersion,  as  we 
shall  see.  Dispersion  refers  to  the  fact  that  surface  waves  of  different  wavelengths  and 
frequencies  travel  at  different  speeds,  and  so  "disperse"  in  space  and  time.  The 
wavelength  and  wave  frequency  are  related  through  a  dispersion  relationship. 

The  dispersion  relationship  is  derivable  from  first  principles:  conservation  of  energy  and 
momentum  for  a  fluid  element  lead  to  the  Navier-Stokes  (NS)  equation  and  the  continuity 
equation.  Detailed  derivations  are  available  in  the  literature  (e.g.  Kinsman  [1965]  or 
Collins  III  [2014]).  Several  assumptions  are  made  about  water  at  this  point  so  that  many 
terms  in  the  NS  are  negligible.  This  includes  neglecting  viscosity  (though  it  is  included  later 
as  one  model  of  attenuation),  compressibility,  and  assuming  water  is  irrotational.  If  water 
is  assumed  to  be  incompressible  and  irrotational,  the  fluid  velocity  can  be  described  by 
the  gradient  of  scalar  function  called  the  velocity  potential.  Without  getting  very  deep 
into  the  details,  the  dynamic  boundary  conditions  at  the  interface  require  the  normal 
stress  to  be  continuous  and  the  shear  stress  to  vanish  (for  inviscid  water  surface).  The 
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stress-strain  relationships  are  determined  by  the  mechanical  model  of  the  ice.  These 
stress-strain  relations  are  then  used  in  the  potential  flow  theory  to  derive  the  dispersion 
relation. 


By  linearizing  the  dynamic  boundary  condition  at  the  surface,  the  open  water  (i.e.  ice  free) 
dispersion  relation  is: 


a)  —  yj  gk  tanh  kd 


For  simplicity,  we  will  take  this  to  the  deep  water  limit,  such  that  tanh  (kd)  goes  to  1: 

.2 


—  or 
=  yjgk  <->  —  =  k  =  kc 
9 


00 


(4) 


(5) 


We  will  maintain  the  deep-water  limit  throughout  the  remainder  of  the  text  though  this 
can  easily  be  extended  to  full  intermediate  water  form.  Here  we  define  kow  as  the  open 
water  wavenumber.  Complications  due  to  ice  are  introduced  through  enforcing  extra 
terms  (i.e.  due  to  the  added  inertia  due  to  the  existence  of  ice  on  the  surface,  the  bending 
due  to  the  flexural  rigidity,  and  the  compressive  stress)  in  the  boundary  conditions  and 
then  re-deriving  the  dispersion  relation.  Details  of  these  derivations  can  be  found 
throughout  the  literature  (e.g.  Squire  et  al.  [1995],  Squire  [2007],  and  Wang  and  Shen 
[2010]). 


2. 1.3.1  MASS-LOADING  MODEL 

The  first,  and  perhaps  simplest,  way  ice  effects  were  introduced  is  by  including  the  added 
mass  at  the  interface;  this  is  known  as  the  mass  loading  model  (originally  developed  by 
Weitz  and  Keller  [1950]). 


g  -  picehoo2/p  (6) 


Following  Liu  and  Mollo-Christensen  [1988],  the  presentation  can  be  simplified  by  defining 
an  inertial  coefficient,  M  —  piceh/p\ 


00 


2 


g  —  Moo2 


—  k 


(7) 


The  model  of  added  mass  is  quite  general,  so  specifics  of  ice  are  introduced  though  the 
inertial  coefficient  with  the  inclusion  of  material  density,  p\ce,  and  thickness,  h.  For 
producing  Figure  2,  the  density  of  sea  water  was  set  as  1025  kgnr3  and  p\ce  was  10%  less 
than  p.  There  is  only  one  free  parameter,  h,  and  as  the  h  goes  to  zero,  the  inertial  term  in 
the  denominator  goes  to  zero  and  the  relationship  simplifies  to  that  of  open,  deep  water 
relationship  of  Eq.  5. 
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Figure  2  The  normalized  dispersion  relation  calculated  from  the  mass  loading  model  as  a  function  of  ice 
thickness,  h.  The  thickness  ranges  from  0.001  to  3.500  meters  linearly  (Ah  =  0.035  m)  where  h  =  0  is  the 
deep  water  linear  dispersion  relation  shown  by  the  dotted  black  line  and  h  =  3.5  meters  is  the  dark  red 
line. 

The  effects  of  different  thickness  ice  on  the  dispersion  relation  is  shown  in  Figure  2.  The 
wavenumber  calculated  with  the  mass  loading  model  is  normalized  by  the  open  water 
wavenumber,  k0w,  so  that  the  y-axis  shows  the  change  in  wavenumber.  Figure  2  shows 
that  mass  loading  model  (1)  always  increases  the  wavenumber  (shortens  the  wavelength) 
and  (2)  that  the  change  in  wavelength  is  most  important  for  higher  frequencies.  Even  for 
very  thin  ice,  h  =  0.035  m,  the  wavenumber  corresponding  to  2  Hz  is  double  that  of  the 
open  water  relationship.  For  long  waves, /<  0.1  Hz,  the  effect  mass  loading  of  ice  up  to 
3.5  m  thick  is  not  substantial.  There  is  a  15%  increase  in  wave  number  at  0.1  Hz. 

The  mass  loading  model  may  be  appropriate  for  conditions  where  the  elastic  properties  of 
ice  are  unimportant  and  floes  do  not  interact  with  other  floes.  This  would  be  the  case  for 
long  waves  propagating  through  small  (compared  to  the  wavelength),  sparse  pancake  ice. 

2. 1.3.2  PURE  ELASTIC  ICE  WITH  MASS  LOADING 

The  next  level  of  complexity  requires  that  some  assumptions  about  the  mechanical 
behavior  of  ice  must  be  made,  which  determine  the  stress  -  strain  relationship.  In  the 
context  of  wave-ice  interaction,  ice  may  be  regarded  as  an  isotropic,  homogenous  elastic 
plate  (or  beam)  due  to  the  strain  rates  in  involved.  The  stress  is  supplied  by  the  incoming 
wave  field,  and  the  strain  describes  the  response  of  the  ice  which,  through  enforcing  the 
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dynamic  boundary  conditions  must  be  continuous  at  the  ice  -  water  interface.  This  is  an 
application  of  the  Euler-Bernoulli  beam  theory  or  its  extension,  Kirchoff-Love  plate  theory. 
Conceptually,  this  might  occur  when  the  ice  is  a  large,  thin  plate  which  surface  waves  may 
induce  coupled  wave  modes.  After  re-deriving  the  dispersion  relation,  the  material 
property  of  flexural  rigidity  shows  up  in  the  denominator  as  a  term  proportional  to  /c4. 

co2 

g  —  Mu)2  +  LkA /p  ^  (8) 

Here  L  is  the  flexural  rigidity  of  the  sea  ice  which  is  a  function  of  three  ice  parameters  (1) 
the  thickness,  h,  (2)  the  effective  elastic  modulus,  E,  and  (3)  the  Poisson  ratio,  v: 

Eh3 

1  ~  12(1  -v2)  (9) 

In  the  purely  elastic  regime,  there  exists  a  simple  relationship  between  the  Lame 
constants  E,  v,  and  G: 

E  =  2G(1  —  v) 

G  =  £7(2(1  -  v)) 


E 


So  we  can  easily  choose  the  free  parameter  to  be  G.  Flexural  term  in  the  dispersion 
relation  tends  to  lengthen  the  wavelength  [Squire,  1993],  The  net  wavelength,  whether 
shorter  or  longer,  depends  on  the  balance  between  mass  loading  and  the  flexural  rigidity. 
We  use  the  Euler-Bernoulli  beam  model  of  Fox  and  Squire  [1994]  (FS)  to  examine  a  few  of 
the  possibilities. 

According  to  Figure  3,  whether  the  wavelength  shortens  or  lengthens  depends  on  the 
relative  contributions  of  the  mass  loading  and  flexural  rigidity  terms  in  the  denominator. 
When  the  shear  modulus  goes  to  zero,  then  the  pure  mass  loading  Eq.  6  is  recovered  from 
Eq.  8. 
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Figure  3  The  normalized  dispersion  relationship  in  ice  with  mass  loading  and  flexural  bending  over  11 
decades  of  shear  modulus.  Ice  thickness  of  0.01  m,  0.1  m,  0.5  m,  1.0  m,  2.0  m,  and  3.5  m  are  shown  from 
left  to  right,  top  to  bottom.  The  shear  modulus  ranges  from  1  Pa  (dark  blue)  to  1010  Pa  (dark  red)  with 
logarithmic  spacing. 

If  we  choose  a  typical  wave  frequency,  say  0.1  Hz,  how  does  the  shear  modulus  change 
the  wavelength?  Figure  4  shows  that  a  0.1  Hz  does  not  change  from  the  mass  loading 
model,  even  in  thick  ice  of  h  =  3.5  m.  For  the  thickness  tested,  deviations  from  the  mass 
loading  model  slowly  appear  beyond  a  shear  modulus  of  106  Pa.  For  a  given  h,  after  the 
effects  of  flexural  appear,  the  wavelength  rapidly  increases  as  a  function  of  G  (and  h).  For 
the  same  level  of  shear  modulus,  say  1010,  a  change  in  ice  thickness  from  1  m  to  2  m  leads 
to  a  large  change  in  wavelength  of  ~60  m. 
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FS  Model  for  0.1  Hz  wave  [h  =  0.01  -  3.5  m] 


5 

o 


G  [Pa] 


Figure  4  Change  in  normalized  wavenumber  as  a  function  of  shear  modulus,  G,  in  the  FS  model.  The  colors 
indicate  ice  thickness  as  indicated  in  the  legend. 

Before  we  move  on  to  dissipative  models  we  should  note  that  the  effects  of  ice 
compression  have  also  been  explored.  These  are  thought  to  be  important  in  pack  ice 
under  high  compressive  stress  [Liu  and  Mollo-Christensen,  1988], 

g  —  Mm2 *  +  Bk4  —  Phk2 /p  ^  (10) 

Where  we  have  introduced  the  bending  coefficient,  B  =  L/p,  and  P  is  the  compressive 
stress.  We  can  further  simplify  by  defining  a  compressive  coefficient,  Q  =  Ph/p2 


2There  is  an  interesting  analogy  between  compressive  stress  in  wave  dispersion  with  the  effects  of  surface 
tension,  the  equation  for  which  bears  the  form: 


&r  =  I gk  + 


r  k'c 


tanh  dk 


Where  r  is  the  surface  tension  of  water.  By  taking  this  form  to  the  deep  water  limit,  and  rearranging: 

0>2 

- rrr  =  k 

g+rk2/p 


Now,  the  surface  tension,  T,  and  the  compressive  stress  can  be  related  by  the  following  equation: 

T 

P  =  ~h 

So  that  compressive  stress  acts  as  negative  surface  tension  scaled  by  the  ice  thickness. 
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(11) 


g  —  Moo2  +  Bk 4  —  Qk2 

The  effect  of  compressive  stress  is  to  increase  the  wavenumber,  but  the  compressive 
stresses  assumed  in  Liu  and  Mollo-Christensen  [1988]  are  considered  unrealistic  under 
most  conditions  [Timco  and  Weeks,  2010]  so  this  formulation  has  lost  favor  in  the 
community.  This  is  perhaps  most  appropriate  during  a  strong  on  ice  storm  event  in  which 
the  ice  edge  becomes  compact  due  to  wave  and  wind  forcing.  Recent  studies  of  ice 
fracture  by  waves  [Asplin  et  al.,  2012;  Kohout  et  al.,  2014;  Collins  et  al.,  2015;  Kohout  et 
al.,  2015]  insinuate  that  perhaps  this  model  deserves  some  renewed  attention. 

2. 1.3. 3  VISCOUS  AND  VISCOELASTIC  MODELS 

Up  until  this  point,  we  have  not  mentioned  wave  energy.  This  is  because  although  the 
previous  models  changed  the  wavelength,  the  wave  energy  is  conserved.  In  contrast, 
introducing  viscosity,  by  itself,  does  not  significantly  change  the  wavelength  for  the 
viscosities  considered  here,  but  it  does  change  the  energy  as  a  function  of  propagation 
distance. 

The  formulation  of  dissipation  may  be  introduced  by  one  of  several  models.  One  may 
introduce  a  viscoelastic  model  such  as  the  Voigt  model  or  Maxwell  model  [Wang  and 
Shen,  2010]  which  determine  the  form  the  deviatoric  stress-strain  relationships.  Different 
physics  are  active  for  different  models  including  a  spring-dashpot  (damped  oscillator)  for 
the  Voigt  model,  viscosity  of  the  water  layer  or  ice  layer  or  both,  viscosity  proportional  to 
the  water  velocity  (friction),  fatigue  or  creep  strain  (one  such  model  is  derived  in  Appeddix 
A:  Bending  and  Dispersion  of  viscoelastic  waves  (Marchenko)),  and  others.  A  class  of  these 
are  known  as  visco-elastic  models  [Wang  and  Shen,  2010;  Mosig  et  al.,  2015], 

To  talk  about  viscous  dissipation  we  need  to  back  up  a  step  and  reformulate  the 
dispersion  in  terms  of  a  complex  wavenumber: 

k  —  k0  +  ia  ^ 

Assuming  the  real  and  imaginary  parts  are  separable3,  solving  for  the  real  part  of  the 
dispersion  relation  gives  the  wavenumber  (as  in  the  previous  examples)  and  solving  for 
the  imaginary  part  gives  the  attenuation  coefficient,  a.  For  a  monochromatic  wave  of 
amplitude  A  which  is  a  function  of  distance,  x: 

A(x )  =  A0elkx 

Substituting  in  the  complex  wavenumber 


3  This  is  apparently  not  possible  for  the  model  of  Wang  and  Shen  [2010]  [J.  Mosig,  personal  communication], 
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A(x)  =  A0eix{-ko+ia)  =  A0e~axeik°x 


Here  the  second  term  on  the  right  hand  side  gives  the  attenuation  and  the  third  term  is 
gives  the  oscillation.  Exponential  attenuation  arises  as  a  natural  consequence  of  the 
complex  wave  number,  this  is  known  in  physics  as  the  evanescent  mode.  Therefore, 
attenuation  is  direct  result  of  solving  the  dispersion  relation.  Given  the  group  velocity, 

Cg  =  the  attenuation  coefficient  may  be  written  in  terms  of  the  space  or  time  domain. 
Giving  the  attenuation  rate  as  a  function  of  length  serves  as  a  convenient  implementation 
of  the  ice  source  term  [Rogers  and  Zieger,  2014]: 


—2  Cga 


(13) 


In  Eq.  11  (and  therefore  its  limits  Eqs.  6  and  8)  dissipation  has  not  yet  been  introduced.  In 
the  non-conservative  regime,  dissipation  follows  directly  from  the  formulation  of  the 
dispersion  relation.  Visco-elastic  models  neglect  compressive  stress  but  expand  on  Eq.  8 
including  an  additional  imaginary  term  within  the  term  proportional  to  k 4. 


oo 


g  —  Moo2  +  DkA 


=  k 


(14) 


Where  D  has  a  real  part  (B  from  above)  and  imaginary  part  consisting  of  any  number  of 
specific  complex  formulations,  one  of  which  we  will  expand  on  below.  We  shall  proceed 
with  an  example  of  a  viscoelastic  model:  the  extended  model  of  Fox  and  Squire,  [1994] 
(EFS).  This  model  follows  from  analogy  to  a  spring-dashpot  and  introduces  a  viscosity  (or 
friction)  proportional  to  frequency  into  the  elastic  plate  model  above.  The  dispersion 
relationship  for  EFS  is  as  follows: 


g  -Moo2  +^f(l  +  v)k4  (15) 


Here  Gv  is  complex  Voigt  shear  modulus  [Mosig  et  al.,  2015]  which  is  related  to  the  elastic 
shear  modulus,  G,  and  the  viscosity  parameter  (related  to  the  dashpot-constant),  4r\\ 

Gv  —  G  i^PiceV  (16) 

So  that  Eq.  15  can  be  expanded  to 


4Not  to  be  confused  with  sea  surface  elevation  for  which  p  is  typically  used  in  the  literature. 

18 


(17) 


co 


g  -  Mo2  +  Bk 4  -  iaopicer]  ^  (1  +  v)k 4 


=  k 


Figure  3  was  produced  by  setting  g  to  zero  in  Eq.  17,  therefore  reducing  the  EFS  to  the  FS, 
and  then  varying  the  value  of  G.  With  all  other  parameters  staying  the  same,  for  the 
values  of  viscosity  tested,  changing  viscosity  alone  does  produce  a  significant  change  to 
the  real  part  of  the  wavenumber  and  hence  does  not  alter  the  wavelength  (this  is  verified 
in  [Mosig  et  al.,  2015]).  It  does  however,  change  the  level  of  attenuation. 


Figure  5  Showing  the  normalized  dispersion  relation  on  the  left  and  the  corresponding  attenuation 
coefficient  on  the  right.  From  top  to  bottom  the  ice  thickness  is  0.1,  0.5, 1.0,  and  3.5  m.  The  colors 
indicate  the  viscosity  ranging  over  5  decades.  The  dispersion  relation  does  not  change  as  a  function  of 
viscosity,  therefore  all  lines  are  over  plotted  and  the  only  visible  line  is  dark  red  (last  one  plotted)  and  the 
black  dotted  line  is  the  open  water  relation. 

Figure  5  shows  the  normalized  dispersion  relation  on  the  left  side  and  on  the  right  the 
attenuation  coefficient.  From  top  to  bottom  the  ice  thickness  is  0.1,  0.5,  1.0,  and  3.5  m, 
respectively.  G  is  set  to  2  x  109to  the  value  used  in  Doble  and  Bidlot  [2013]  and  the 
variation  of  v  is  shown  with  color.  Because  dispersion  is  not  a  function  of  viscosity  in  this 
formulation,  all  the  dispersion  curves  lie  on  top  of  one  another,  leaving  only  the  most  last 
one  plotted  (dark  red)  visible.  In  contrast,  the  attenuation  rate  varies  monotonically  with 
viscosity.  The  change  in  slope  in  attenuation  space  can  be  seen  to  correspond  with 
deviation  from  the  open  water  dispersion  relation.  Attenuation  rates  are  highest  for  the 
high  frequencies.  Given  the  same  viscosity,  whether  or  not  attenuation  increases  or 
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decreases  with  a  change  in  ice  thickness  depends  on  the  frequency.  Comparing  the  top 
and  bottom  plots  on  the  right  hand  side,  the  lower  frequency  (<  0.10  Hz)  attenuation  is 
drastically  increased  and  the  higher  frequency  (>  0.30  Hz)  attenuation  is  slightly 
decreased. 

To  explore  this  a  bit  further,  we  set  reasonable  values  for  elastic  shear  modulus,  G  =  105, 
and  viscous  parameter,  77  =  0.01.  We  vary  the  ice  thickness  from  0  -  5  m,  the  value  of 
thickness  indicated  by  color  ranging  from  dark  blue  to  red,  respectively.  For  reference, 
solution  for  the  dispersion  relation  with  ice  thickness  of  5  m  for  the  mass  loading  model 
(ML)  and  purely  elastic  model  (FS). 


G  =  105  [Pa];  v  =  0.01  [mV1];  H  =  0.0  -  5.0  A0.5  [m] 


Figure  6  Left  side:  Dispersion  relation  according  to  the  extended  Fox  and  Squire  model.  The  color  indicates 
ice  thickness,  h.  Shear  modulus  is  set  to  the  value  10s  and  viscous  parameter  to  the  value  0.01.  The  mass 
loading  model  (ML)  and  purely  elastic  model  (FS)  are  shown  with  the  black  dotted  line  and  black  dashed 
lined,  respectively.  Right  side:  Attenuation  coefficient  corresponds  to  the  left  side. 

The  left  of  Figure  6  shows  the  dispersion  relation  calculated  using  the  EFS  model  with 
color  showing  ice  thickness.  The  solution  for  h  =  5,  matches  that  of  the  purely  elastic 
model  (FS)  exactly.  All  three  models  (EFS,  FS,  and  ML)  give  indistinguishable  k  for 
frequencies  less  than  0.15.  On  right  side,  attenuation  is  a  function  of  frequency,  with 
attenuation  most  important  for  the  high  frequencies.  The  transition  in  slope  for  the 
attenuation  corresponds  with  the  dispersion  relation  transitioning  from  mass  loading 
dominant  to  elasticity  dominant.  Before  the  transition,  the  attenuation  increased 
monotonically  with  thickness,  the  opposite  is  true  after  the  transition,  resulting  in  an 
unintuitive  situation:  for  high  frequency  (>  0.5  Hz)  waves  less  ice  means  more  attenuation. 
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That  being  said,  one  should  keep  in  mind  that  relationship  being  tested  is  based  on  Euler- 
Bernoulli  beam  theory  which  valid  only  for  infinitesimal  strains.  For  high  frequency  waves 
(>  0.5-  Hz),  the  reflection  will  be  nearly  100%  [Fox  and  Squire,  1995]  so  a  phenomenal 
amplitude  is  required  to  register  strain.  In  essence,  such  a  situation  is  not  physically 
feasible  to  begin  with. 

Lastly  there  are  some  artifacts  to  point  out.  Where  a  dispersion  relation  for  a  particular 
thickness  appears  to  drop-off  or  disappear  is  an  error  (correspondingly  in  attenuation 
appears  step-wise).  The  solver  is  giving  the  wrong  root  (e.g.  see  Mosig  et  al.  [2015])  and 
the  parameters  involved  are  outside  the  validity  of  the  model  (personal  communication 
with  V.  Squire  [2015]).  Apparently,  compared  to  the  model  of  Wang  and  Shen  [2010],  the 
EFS  model  is  the  less  susceptible  to  these  sorts  of  errors  and  great  care  needs  to  be  taken 
to  ensure  the  proper  root  is  chosen  (again,  see  Mosig  et  al.  [2015]).  Devising  a  system  for 
choosing  roots  of  the  complex  dispersion  relation  (e.g.  Zhao  et  al.  [2015])  is  an  active  area 
of  research. 

In  the  limit  that  ice  is  homogeneous,  one  could  specify  which  dispersion  relation  is 
appropriate  depending  on  the  ice  and  wave  properties.  So  the  running  assumption  is  an 
isotropic  ice  medium.  Anisotropies  and  organized  ice  patterns  will  lead  to  deviations  from 
the  derivations  below.  For  large-scale,  but  not  fine-scale  homogeneity,  different 
properties  of  the  ice  will  manifest  in  the  measured  dispersion  relation.  Using  an  average 
ice  thickness  where  in  reality  ice  thickness  varies  would  be  analogous  estimating  a 
vehicle's  speed  by  dividing  the  distance  covered  by  the  travel  time,  when  in  reality,  for 
whatever  reason  (perhaps  traffic),  the  speed  of  the  vehicle  was  not  constant.  In  this  way, 
one  might  refer  to  a  measured  dispersion  relation  as  a  bulk  dispersion  relation.  This  being 
said,  the  assumption  of  ice  as  a  continuum  may  be  a  gross  oversimplification  in  some 
cases. 

2. 1.3.4  DISPERIONS  MEASUREMENTS 

We  have  introduced  several  simple  mathematical  treatments  for  dispersion  of  surface 
waves  in  ice,  and  although  equations  are  somewhat  simple  to  show,  the  solutions 
produced  are  very  mathematically  rich.  The  measurements  of  the  dispersion  relationship 
in  ice  are  much  fewer  than  those  of  attenuation,  and  hence  bulk  dispersion  relationship  is 
even  less  well  understood  than  bulk  attenuation  in  practice.  The  difficulty  is  that 
estimation  of  the  dispersion  requires  spatio-temporal  information,  e.g.  co-located 
measurements  of  wave  period  and  wave  length.  In  the  end,  any  measurements  of 
dispersion  in  ice  are  a  manifestation  of  many  individual  mechanisms  of  wave  ice 
interaction  and  thus,  like  attenuation,  is  bulk  property.  Whereas  the  mechanisms 
described  above  are  predominately  dissipative  and  clearly  lead  to  bulk  attenuation,  the 
dispersive  outcomes  of  the  individual  mechanisms  are  mixed  (could  be  a  balance  of 
mechanisms  which  increased  and  decreased  wavelength)  leading  to  ambiguous  results. 
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Several  studies  attempted  to  use  SAR  data  to  measure  the  dispersion  relation  [Wadhams 
and  Holt,  1991;  Liu  et  al.,  1991],  For  the  most  part,  it  was  found  that  the  wavelength 
decreased  which  supported  the  mass  loading  model  in  the  case  of  Wadhams  and  Holt 
[1991]  and  high  compression  in  the  case  of  Liu  et  al.  [1991],  though  evidence  was 
generally  limited  in  both  studies.  Fox  and  Haskell  [2001]  was  able  to  estimate  the 
propagation  speeds  of  waves  (and  hence  indirectly  the  wavelength)  in  ice  by  measuring 
the  frequency  spectrum  at  two  closely  located  positions.  Their  Figure  6  shows  that  fitted 
empirical  wavelength  is  slightly  longer  for  frequencies  within  the  0.05  -  0.10  Hz  band  and 
then  significantly  shorter  for  frequencies  from  the  0.10  -  0.16  Hz  band.  Sakai  and  Hanai 
[2002]  in  a  lab  experiment  with  synthetic  ice,  found  that  the  dispersion  relationship 
(between  0.6  and  1.7  Hz)  varied  between  a  flexural  model  and  a  mass  loading  model  as  a 
function  of  floe  length  scale.  The  transition  must  have  been  dependent  on  the  relative 
scales  of  the  wavelengths  and  floe  sizes  so  that  at  any  one  time  the  measured  dispersion 
relation  is  a  mixture  of  the  mass  loading  and  flexural-gravity  models.  The  higher 
frequencies  showed  the  greatest  celerity  ranges  (the  range  itself  a  function  of  ice 
thickness). 

Under  very  precise  laboratory  conditions,  literally  performed  on  a  desk  top,  both  the 
dispersion  relationship  and  attenuation  has  been  directly  measured  in  viscous  water  for 
waves  in  the  gravity-capillary  wave  range  [Behroozi  et  al.,  2010],  By  controlling  the  wave 
frequency  and  generating  a  standing  wave,  they  were  able  to  measure  the  wavelength 
(giving  dispersion)  and  the  change  in  wave  height  with  space  (attenuation).  They  were 
then  able  to  invert  the  dispersion  relation  to  give  an  estimate  of  the  viscosity.  In  the 
context  of  field  experiments,  there  are  obvious  technical  challenges  which  differentiating 
the  two  settings:  a  spectrum,  not  a  single  frequency,  of  waves  will  be  freely  propagating, 
not  standing.  A  standing  wave  removes  the  spatial  dependence  of  wave  properties,  so  in 
the  field  a  spatio-temporal  measurement  needed,  of  which  there  are  very  few  examples 
and  none  in  ice. 

Here  we  speculate  on  how  to  accomplish  such  a  measurement.  There  are  few  techniques 
for  spatio-temporal  field  measurements,  one  category  is  stereo-video  systems  (e.g. 
Campbell  et  al.  [2014])  and  another  is  ship  borne  X-band  marine  radar  systems  (MR).  We 
concentrate  on  the  latter  since  MR  has  recently  been  utilized  in  icy  conditions  as  part  of 
the  Sea  State  DRI.  MRs  transmit  microwave  pulses  and  record  the  backscatter  intensity, 
and  backscatter  is  dominated  by  the  Bragg  scattering  mechanism  [Young  et  al.,  1985; 
Borge  et  al.,  1999],  The  backscatter  is  tracked  in  space  and  time.  The  time  of  return 
provides  the  radial  distance  (range)  and  the  look  direction  of  the  antenna  provides  the 
angle.  The  rate  of  rotation  can  be  such  that  MRs  produce  one  complete  (360°)  scan  of  the 
sea  surface  nearly  every  second.  One  then  performs  a  3D  Fourier  transform  of  the  time 
series  snapshots  to  get  a  3D  image  spectrum  and  which  is  then  converted  to  a  2D 
wavenumber  spectrum,  details  for  which  may  be  found  in  (e.g.  Lund  et  al.  [2014]). 
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The  cornerstone  to  the  MR  method  is  a  direct  measurement  of  the  location  of  wave 
energy  in  3D  Fourier  space  which  is  defined  by  the  dispersion  relation: 

co  =  F{k )  +  k  U 

Where  F(k)  is  whatever  form  of  the  dispersion  relationship  is  applicable  and  U  is  the 
surface  current  vector.  Usually,  the  dispersion  relation  for  linear,  open  water  is  the  only 
applicable  dispersion  relation,  so  that  deviations  from  the  linear  dispersion  curve  indicate 
the  presence  of  a  current.  A  wave  spectrum  under  the  influence  of  a  uniform  surface 
current  experiences  a  Doppler  shift.  Therefore,  measured  deviations  from  the  linear 
dispersion  relation  are  assumed  to  be  in  the  form  of  a  current-induced  Doppler  shift  (e.g. 
Senet  et  al.  [2001]). 


Figure  7  Reproduced  from  Figure  4  A.  of  Lund  et  al.  [2015]  (copyright  American  Meteorological  Society) 
see  the  reference  for  full  details.  Cross-section  through  the  3D  image  spectrum  in  the  peak  wave 
direction.  Wavenumber  on  the  x-axis  and  angular  frequency  on  the  y-axis.  "Still  water"  represents  the 
open  water,  linear  dispersion  relation  (without  current)  shown  by  a  solid  black  line.  The  color  code  (blue 
to  red)  indicates  the  energy  detected  by  MR.  A  speed  of  0.52  ms  1  produces  a  good  fit  (black  dashed  line) 
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to  the  dispersion  curve  measured  by  the  radar  (energy  significantly  above  the  background  noise  level). 
Higher  order  harmonics  are  detectable  in  the  MR  spectrum. 

Figure  6,  reproduced  from  Lund  et  al.  [2015],  shows  an  example  of  this  process.  It 
represents  a  cross-section  through  one  3D  image  spectrum  in  the  peak  wave  direction. 
Here  we  have  wavenumber  on  the  x-axis  and  angular  frequency  on  the  y-axis,  we  will 
refer  to  this  as  the  dispersion  plane.  The  colors  in  the  background  correspond  to  energy 
mapped  to  dispersion  plane.  The  solid  black  line  labeled  "Still  water"  is  the  open  water, 
linear  dispersion  relation  (without  current).  Clearly,  the  energy  detected  by  MR  does  not 
follow  this  curve.  For  a  particular  angular  frequency,  the  wave  numbers  are  lower  (the 
wavelength  is  increased),  from  what  one  would  expect  from  the  linear  dispersion  relation. 
Using  an  iterative  fitting  process,  the  surface  current  which  would  result  in  such  a  Doppler 
shift  matching  the  location  of  the  energy  is  found.  In  this  way,  the  surface  current  is 
always  a  byproduct  of  MR  wave  analysis.  For  the  example  above,  a  current  of  0.52  ms  1 
provided  a  Doppler  shifted  dispersion  curve  which  was  a  good  fit  to  the  location  of  the 
energy  on  the  dispersion  plane. 

In  ice  the  dispersion  relation  will  become  the  unknown,  so  an  independent  current 
measurement  will  be  necessary.  With  U  known,  one  can  predict  the  expected  Doppler 
term  and  hence  deviation  from  the  linear,  open  water  dispersion  relation  so  that  any 
further  deviation  must  be  due  to  ice  cover.  Even  in  the  absence  of  an  independent  current 
measurement,  one  might  be  able  to  make  a  rough  guess  at  the  dispersion  relation  if  the 
currents  calculated  by  MR  are  unreasonable. 

Continuing  the  discussion  on  MR  measurement,  once  the  wave  spectra  are  calibrated,  MR 
may  provide  a  spatio-temporal  approach  for  calculating  attenuation.  The  former 
technique  requires  open  water  and  ice  in  a  single  image.  If  this  is  the  case  one  could  use 
different  analysis  boxes  on  a  single  image,  marching  into  the  ice,  to  determine 
attenuation.  The  main  problem  with  this  approach  is  that  there  is  a  known  dependency  of 
wave  results  on  MR  range  and  look  direction  [Lund  et  al.,  2014],  A  more  robust  approach 
would  be  to  use  the  optimal  field  of  view  and  over  time  go  from  open  water  to  icy 
conditions. 

Once  dispersion  and  attenuation  are  known,  and  if5  the  other  measurable  properties  of 
ice  in  visco-elastic  models  (ice  thickness  and  elastic  modulus)  are  known,  then  the  visco¬ 
elastic  dispersion  relation  can  be  inverted  to  give  effective  viscosity  which  is  otherwise 
unmeasurable.  This  would  be  quite  a  feat. 


5This  is  a  huge  if  as  there  a  many  uncertainties  related  to  the  mechanical  properties  of  ice  (see  Timco  and 
Weeks  [2010]). 
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2.1.4  Refraction  and  Shoaling 

One  of  the  important  consequences  of  an  altered  dispersion  relationship  is  the  possible 
change  in  wave  height.  In  direct  analogy  for  refraction  and  shoaling  in  shallow  water  (or  in 
currents),  it  can  be  shown  that  a  change  in  group  speed  necessarily  requires  a  change  in 
direction  and  amplitude  (this  is  Snell's  law,  most  commonly  encountered  in  the  context  of 
geometric  optics).  Furthermore,  if  the  wavelength  shortens,  then  the  waves  will  turn, 
similar  to  a  situation  of  waves  approaching  the  shore:  waves  always  turn  to  approach 
shore-normal  (and  vice  versa  if  the  wavelength  elongates  in  ice).  Let's  imagine  a  situation 
where  monochromatic  waves  of  wavelength  A  approach  an  ice  field  with  properties  which 
slowly  vary  over  uniform  contours.  Wave  crests  approach  the  ice  field  some  angle,  0,  off- 
normal  to  the  ice  contour  parallels  and  let  4  be  the  corresponding  angle  between  the 
wave  crests  and  contour  parallels.  According  to  Snell's  law,  waves  traveling  in  material  1 
and  2  must  obey  the  following  relationship: 


sin  (p± 
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Where  c  is  the  group  speed.  The  group  speed  in  open,  deep  water  (ow)  follows  from  the 
linear  dispersion  relation: 
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As  waves  move  into  ice  covered  regions,  Snell's  law  will  determine  the  refraction  angle: 

sin  iCe  sin  @ice  Cgice  kow  Lice 
sin  (pow  sin  90w  CgOW  k[ce  Low 

The  change  in  wave  height  is  derived  by  conserving  the  mean  energy  per  unit  area,  which 
can  be  determined  from  the  spectrum  by  linear  theory  [Kinsman,  1965]: 


E  —  2pg  j  S(f)df  =  ^pgH^0 


The  wave  energy  is  simply  the  spectrum  scaled  by  gravity  and  the  density  of  water  and  the 
wave  energy  propagates  at  the  group  speed 

P  —  ECn 


Assuming  that  the  energy  contained  between  two  parallel  rays  is  conserved  along  a 
section  of  crest  width,  s. 
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Where  D  is  the  square  root  of  ratio  of  the  group  velocities,  which  the  appropriate  cg,iCe 
may  be  obtained  from  above. 
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K  is  the  square  root  of  the  ratio  of  the  crest  lengths,  from  Snell's  law  above 

i 
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The  implication  is  this:  if  the  group  velocity  slows  in  ice,  then  wave  heights  will  increase 
(shoal)  and  if  the  group  velocity  increases  in  ice,  then  wave  heights  will  decrease.  In  terms 
of  the  wave  effects  on  ice,  this  will  change  the  strain  felt  by  the  ice.  Let  say  the  group 
velocity  decrease  increasing  the  wave  height:  because  of  the  increased  wave  height,  ice- 
shoaled  waves  will  put  a  greater  (than  expected  from  open  water)  strain  on  the  ice 
leading  to  breaking  events  which  might  not  otherwise  occur.  Over  long  distances, 
attenuation  must  be  considered,  and  the  net  wave  height  will  be  a  result  of  shoaling 
balanced  against  energy  lost  to  dissipation  or  scattering.  In  addition,  shortened 
wavelength  and  increased  wave  height  increases  steepness  of  the  waves.  The  steepness 
may  increase  the  likelihood  that  the  waves  will  be  unstable  to  modulational  perturbations 
or  breaking  (see  section  4.1),  this  has  already  been  shown  to  be  important  for  waves 
shoaling  on  an  opposing  current  [Toffoli  et  al.,  2015], 


2.2  WAVE  EFFECTS  ON  ICE 


2.2.1  Overview 

This  document  is  admittedly  biased  towards  ice  effects  on  waves  reflecting  the  first 
author's  background  in  waves.  There  is  rich  body  of  research  on  ice  mechanics,  of  which 
the  effect  of  waves  is  only  one  of  many  possible  influences.  Here  we  attempt  to 
summarize  the  more  salient  features  of  wave  effects  on  ice. 
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Many  wave  effects  on  ice  were  previously  mentioned.  A  short  list  includes  the  response  to 
wave  induced  stress:  recoverable  strain  (elasticity),  elastic  hysteresis,  deformation 
(plasticity),  creep,  creep  buckling  (failure),  elastic  buckling  (failure),  creep  crushing, 
continuous  crushing,  horizontal  splitting,  radial  cracking,  bending  failure,  non- 
simultaneous  failures  [Frannson,  2009],  Although  this  is  a  long  list  of  specific  terms,  the 
connection  between  the  wave  field  and  these  mechanics  is  often  tenuous. 

The  most  important,  in  terms  of  melting,  is  facture  of  the  ice.  Fracturing  results  in  a 
geometrical  or  topological  change  which  effectively  increases  the  lateral  surface  area  of 
the  ice  and  accelerates  the  melting  process  through  this  increased  exposure  [Squire, 

2007],  A  much  cited  reference  for  lateral  melting  is  Steele  [1992],  who  used  a  simple  ice- 
ocean  model  to  show  that,  for  conditions  similar  to  those  of  summer  time  ice  retreat, 
lateral  melt  was  significant  only  for  floes  with  diameters  less  than  0(30  m).  We  should 
note  here  that  the  thermodynamics  in  this  model  are  based  on  diffusion,  that  heat  flux  is 
proportional  to  surface  area  and  that  larger  scale  ocean  dynamics,  such  as  eddies,  are  not 
accounted  for.  Though  the  largest  role  played  by  the  waves  is  increasing  lateral  surface 
area  through  breaking,  it  is  also  mentioned  that  waves  in  leads  "can  accelerate  lateral 
melting  near  the  surface"  presumably  due  to  mixing.  Additionally,  waves  may  over  top 
and  wash-over  ice  surfaces,  which  is  currently  being  investigated  as  a  new  mechanism 
contributing  to  ice  melting  and  reduced  wave  transmission  [Skene  et  al.,  2015;  Toffoli  et 
al.,  2015], 

2.2.2  Modeling  Ice 

Over  an  area  of  sea  surface,  the  ice  can  be  described  by  one  or  more  of  several 
properties:  ice  thickness  distribution  (ITD),  floe  size  distribution  (FSD),  and  ice 
concentration  (1C).  The  ITD  describes  the  percent  surface  area  of  ice  attributable  to  a 
particular  ice  thickness,  h.  Similarly  the  FSD  describes  the  percentage  of  surface  area  of 
ice  attributable  to  ice  of  a  diameter  (some  treatments  consider  caliper  length),  /.  The 
concentration  is  the  percent  area  covered  by  ice  and  typically  acts  to  scale  the  influence 
of  the  ice. 

CICE  (Los  Alamos  Community  Sea  Ice  Model)  [Hunke  et  al.,  2010]  is  a  sophisticated, 
operationally  run  ice  model  based  upon  the  continuity  equation  for  ITD  first  described  in 
Thorndike  et  al.  [1975],  Waves  have  not  played  a  role  in  CICE  because  their  influence  is 
limited  to  the  MIZ  which  was  traditionally  a  small  part  of  the  bulk  ice  field.  However,  the 
MIZs  are  increasingly  accounting  for  higher  percentages  of  the  total  ice  coverage 
particularly  in  the  summer  months,  and  now  account  for  >50%  [Strong  and  Rigor,  2013], 
How  to  get  wave  information  into  CICE,  and  what  to  do  with  it  once  it  is  there,  is  a  topic  of 
active  research.  By  contrast,  prognostic  models  based  on  continuity  of  FSD  and  joint  floe 
size  and  thickness  distribution  (FSTD)  are  a  very  new  development  but  wave  physics  arises 
naturally  as  a  key  player  driving  the  their  evolution  [Zhang  et  al.,  2015;  Horvat  and 
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Tziperman,  2015],  One  path  forward  might  be  to  use  a  FSD  model,  or  breaking  algorithm, 
as  a  mediator  between  a  wave  model  and  an  ITD  model. 


We  now  dive  a  little  deeper  into  the  inner  workings  of  an  ice  model  by  deriving  a 
continuity  equation  using  ice  property,  p,  which  could  be  caliper  length,  diameter  (for 
round  floes),  thickness,  or  some  combination  of  these.  The  generalized  distribution  (which 
could  be  ITD,  FSD,  or  FSTD)  is  then  d( p).  Following  Thorndike  et  al,  [1979]  d(p)  is  area 
conserving,  so  that  the  integral  over  all  ranges  of  ice  property,  p,  is  unity. 

rPmax 

|  d(p)dp  —  1 

J  o 


And  the  area.  A,  covered  by  ice  within  the  property  range  px  <  p  <  p2 relative  to  the 
total  area,  R,  is 

rP2  1 

d(p)dp  =  -A(p1,p2) 


The  distribution  evolves  in  space  and  time  according  to  the  prognostic  continuity  equation 
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This  is  the  familiar  fluid  continuity  equation,  where  the  left-hand  side  is  the  change  in  time 
and  the  first  term  on  the  right-hand  side  is  advection.  In  the  second  term  on  the  right- 
hand  side,  L  represents  the  different  sources,  sinks,  and  redistributors  of  the  distribution 
of  ice  property.  L  consists  of  terms  which  parameterize  the  effects  of  mechanical 
interaction,  thermodynamics,  ridging/rafting  and  among  others.  Very  recently,  some  work 
has  been  done  to  add  a  term  to  L  which  parameterizes  the  wave  effects  on  FSD  and  FSTD 
which  we  will  take  up  in  the  following  sections. 


For  completeness  we  note  that  d(p),  e.g.  FSD,  is  not  consistently  defined  in  the  literature. 
Flere  we  have  based  the  definition  on  area  distribution  in  analogy  with  ITD.  Also  in  use  is 
number  distribution,  N(p),  and  cumulative  number  distribution  where 


N(p)  = 


d(_p ) 

np2 


2.2.3  Floe  Size  Distribution 

Floe  size  distribution  describes  the  number  of  floes  of  different  sizes  over  an  area.  FSD  is 
an  important  concept  in  the  MIZ,  as  the  pack  ice  in  the  central  area  of  the  Arctic  is  fairly 
uniform  with  large,  indistinguishable  floes  with  little  or  no  open  water.  In  the  MIZ,  as  one 
goes  from  open  water  towards  the  pack  ice,  the  FSD  is  characterized  by  large  number  of 
small  floes  with  more  open  water  and  then  a  smaller  number  of  larger  floes  with  less 
open  water. 
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The  exceedance  distribution  (probability  that  a  floe  exceeds  a  particular  size)  is  defined  as 

P(Z  >  Z*)  =  1  -  f  d(p)dp  for  l  >  lmin 
Jo 

Here  /  is  a  parameter  which  describes  the  size  of  the  floe  and  /min  is  the  minimum  floe  size 
considered.  In  reality,  there  is  no  limit  to  the  minimum  floe  size,  but  practically,  the 
minimum  observable  flow  size  depends  on  the  resolution  of  the  image.  And  the  minimum 
floe  size  considered  in  primitive  wave-ice  interaction  models  is  typically  0(20  m)  so  that 
scattering  based  wave  attenuation  (which  increases  strength  with  the  number  of  floes) 
doesn't  become  overly  strong.  This  is  obviously  a  problem  for  a  coupled  thermodynamic 
model  because  this  leaves  little  overlap  between  minimum  modeled  floe  size  0(20  m)  and 
the  floe  sizes  most  effected  by  lateral  melt  <  0(30  m). 

Observations  of  FSD  [Rothrock  and  Thorndike,  1984;  Toyota  et  al.,  2006],  done  through 
analysis  of  aerial  images,  have  shown  exceedance  distribution  follows  a  power-law 
distribution  with  exponent,  —  y.  A  power-law  distribution  implies  scale  invariance  or 
fractal  nature.  Several  studies  have  shown  y  ~  2  [Rothrock  and  Thorndike,  1984;  Perovich 
and  Jones,  2014],  Toyota  et  al.  [2006]  observed  y  —  1.87  (also  gives  a  range  from  1  to  5) 
for  floes  greater  than  0(40  m)  and  y  =  1.15  for  floes  less  than  0(40  m).  Apparently,  wave 
induced  breaking  was  the  dominant  mechanism  for  setting  the  exponent  for  smaller  floes. 
Perovich  and  Jones  [2014]  observed  an  altogether  departure  from  the  power  law  behavior 
for  smaller  floes,  hypothesized  to  be  due  to  lateral  melting.  The  consensus  on  a  power  law 
distribution  is  tenuous:  as  shown  by  Herman  [2010],  the  scatter  from  one  dataset  to 
another  is  considerable  and  may  be  better  fit  by  other  distributions. 

2.2.4  Primitive  Wave-Ice  Interaction  Models 

At  the  time  of  writing,  there  are  only  a  handful  of  FSD  models  which  handle  coupled 
wave-ice  interaction,  so  that  we  try  to  summarize  the  approach  of  each.  There  are  two 
schools  of  thought,  one  which  aims  to  develop  a  redistribution  term  for  prognostic 
prediction  of  the  FSD,  and  another  which  takes  advantage  of  the  parametric  power-law 
description  of  FSD.  Starting  with  the  former,  Thorndike  et  al.  [1975]  first  used  the  general 
equation  for  non-local  interaction  of  a  property  distribution  for  ITD  due  to  ridging 

4  =  —Q(p)d(p)  +  I  (3(p'  ,p)Q(p')d(p')dp' 

Jo 

Zhang  [2015]  followed  this  approach  for  a  redistribution  function  for  FSD  due  to  stochastic 
wave  forcing.  The  ice  property  is  caliper  length.  /?  is  the  redistributor  of  FSD  and  Q  is  the 
redistribution  probability  function.  Under  the  constraint  that  /°°  /?(p' ,p)dp  =  1, 

n,  ,  =  fl/(C2Pl  -  CiPi),  if  CiPi  <  p2  <  C2p1 

PK.P  ,PJ  [0,  i/  l2  <  c-lP-l  or  p2  >  C2P-L 
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ci  and  C2  are  redistribution  coefficient  constants  defined  by  the  minimum  and  maximum 
caliper  size.  The  redistribution  probability  function  is  as  follows: 


Q(p )  =  max 


0 


The  open  parameter  here  is  Cb,  the  participation  factor,  which  in  general  takes  on  a  value 
1  >  cb  >  0.  This  set  up  ensures  that  when  floe  sizes  break,  there  is  equal  probability  that 
any  floe  size  bin  smaller  than  the  original  floe  bin  size  may  gain  from  breaking  (see  the  full 
reference  Zhang  [2015]  for  details).  The  result  of  breaking  is  then  a  FSD  which  obeys  a 
power  law  distribution  and  puts  the  model  in  qualitative  agreement  with  this  feature  of 
FSD  observations. 


Another  approach  is  the  wave  focused  studies  of  Williams  et  al.  [2013a,  b]  which  is  based 
on  the  work  of  Dumont  et  al.  [2011],  Flere,  a  power  law  FSD  is  assumed,  breaking  then 
only  adjusts  the  mean  and  maximum  floe  sizes.  The  minimum  floe  size  is  set  to  20  m  and 
the  mean  floe  size  is 


(Dn) 


y 

y  —  n 


Dn 

^min 


(18) 


Flere  y  is  related  to  a  breaking  probability  FI  that  a  floe  will  break  into  <f2  pieces 

y  =  2+  logf  77 

11  is  set  to  0.9  and  ^  to  2.  This  relationship  and  the  values  chosen  are  based  on  previous 
work  by  [Toyota  et  al.,  2011],  This  effectively  gives  y  —  1.848. 

The  model  was  motivated  by  the  following  scenario:  imagine  waves,  as  predicted  by  a  3G 
wave  model,  approaching  the  MIZ.  Before  entering  the  icy  grid  point,  this  model  would 
take  over  and  predict  1)  the  wave  spectrum  in  icy  waters  and  2)  the  FSD  throughout  the 
MIZ  subject  to  waves.  Given  an  initial  FSD  and  initial  wave  spectrum,  the  spectrum  is 
marched  into  icy  grid  point  where  it  is  subject  to  attenuation  due  to  the  initial  ice 
conditions.  The  ice-effected  wave  spectrum  is  used  to  calculate  a  significant  strain  on  the 
ice  based  on  plate  theory: 


Flere  h  is  ice  thickness  and  qICe  is  the  ice  surface  elevation.  Since  a  wave  spectrum  is 
stochastic  (random  phase)  and  not  deterministic,  the  strain  exceedance  probability 
function  is  developed  based  on  the  Rayleigh  distribution: 

Fe  =  P(e  >  £c)  =  exp(—£c /2m0[s]) 

Where 
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conS(a)) 


h 


.2 
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dw 


3(<w)  is  a  transmission  coefficient  and  the  critical  strain  is  assumed,  ec  —  3  x  10-5.  Given 
the  length  of  time  step  and  the  zero-crossing  wave-in-ice  period,  then  the  number  zero- 
crossings  is  known  and  the  exceedance  probability  may  be  evaluated.  If  critical  strain  level 
is  exceeded,  the  floes  break  and  then  maximum  floe  size  is  set  to  half  the  dominant 
wavelength.  Squire  and  Montiel  [2015]  offer  small  improvements  by  giving  a  nuanced 
definition  for  ice  strength  and  including  vertical  variations  of  ice  properties. 


Horvat  and  Tziperman  [2015]  developed  a  prognostic  equation  for  joint  FSTD,  similar  in 
vein  to  Zhang  [2015]  but  with  explicit  physics  for  the  processes  of  freezing/melting, 
ridging/rafting,  and  breaking  due  to  waves.  The  treatment  of  waves  follows  Williams  et  al. 
[2013a],  but  lacks  some  of  the  nuanced  details.  For  example,  instead  of  using  a  visco¬ 
elastic  dispersion  relation  to  derive  wavelength,  the  open  water  dispersion  relation  is 
used.  The  model  is  quite  sophisticated  in  other  ways,  and  summarizing  the  rest  of  the 
model  is  outside  the  scope  of  this  study  but  we  encourage  the  interested  reader  to  see 
the  full  paper. 


To  summarize,  the  following  table  describes  the  treatment  of  waves  and  ice  published 
wave-ice  interaction  models. 
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Table  1  Summary  of  primitive  wave-ice  interaction  and  FSD  models. 


Study 

Wave  Treatment 

Floe  Treatment 

Strain/Wave 
Height  from 
Spectrum 

Dispersion 

Attenuation 

Fracture 

FSD 

calculation 

Kohout  and 
Meylan  [2008] 

Deterministic 

Elastic  (e.g.  Fox  and 
Squire  [1994]) 

Scattering 

Stress  Criteria 

N/A 

Dumont  et  al. 
[2011] 

Deterministic 

Open  Water 

Scattering  (Kohout 
and  Meylan  2008) 

Strain  and  Stress 

Criteria 

Constrained  by 
power  law  model 
(Eq.  18) 

Williams  et  al. 
[2013a,  b] 

Probability 

Distribution 

Function 

Visco-elastic 
(Robinson  and 
Palmer  [1990]) 

Viscous  (Robinson 
and  Palmer  [1990]) 
and  Scattering 
(Bennets  and 
Squire,  [2012]) 
(scattering 
dominant) 

Stress  Criteria 

Constrained  by 
power  law  model 
(Eq.  18) 

Zhang  [2015] 

N/A 

N/A 

N/A 

Parameterized 

Constant 

distribution  with 
parameterized 
min/max 

Horvat  and 
Tziperman 
[2015] 

Stochastic 

Open  Water 

Parameterized 
Scattering  based  on 
a  fit  to  Kohout  and 
Meylan  [2008] 

Stress  Criteria 

From  wave 

information 

Clearly  there  is  much  room  for  improvement  and  experimentation.  On  the  side  of  wave 
physics,  wave  height  is  a  stochastic  process,  so  the  probability  of  ice  breaking  is  related  to 
the  wave  height  exceedance  probability,  in  this  case  the  Rayleigh  distribution  has  been 
used  extensively.  The  Rayleigh  distribution  was  derived  for  linear,  narrow  banded  wave 
spectrum.  A  simple  improvement  of  the  model  would  be  to  some  other  distribution  (e.g. 
Tayfun  and  Fedele  [2007]),  but  this  is  a  fine  tuning  compared  to  the  gross  uncertainty  in 
ice  properties.  For  these  models  to  significantly  improve,  we  need  much  better 
observations  of  the  mechanical  properties  of  ice  especially  in  concert  with  detailed  wave 
measurements.  We  briefly  address  ice  in  the  following  section. 

2.2.5  Ice  Response 

If  the  response  of  ice  to  stress  is  purely  elastic  (linear)  a  generalization  of  Flooke's  law  may 
be  applied.  Flere  the  stress  and  strain  are  proportional  with  a  coefficient  known  as  the 
effective  elastic  modulus  (Young's  modulus).  Under  these  linear  conditions,  the 
relationship  between  Lame  constants  Young's  modulus,  the  strain  modulus,  Poisson's 
ratio  is  derivable: 


E  =  2G(1  —  v) 

G  —  £7(2(1  —  v)) 


Typically,  the  strain  response  of  ice  is  not  purely  elastic  either  due  to  stress  level  effects 
such  as  yield  or  due  cumulative  stress  effects  such  as  creep  strain  [Timco  and  Weeks, 
2010],  Under  these  conditions,  the  relationships  above  are  not  so  clear,  and  the  modulus 
relationships  no  longer  describe  the  strictly  linear  portion  of  the  stress  strain  relationship. 
Stress  -  strain  relationship  may  still  be  described  with  effective  modulus  which  include 
strain  due  to  creep  and  inelasticity.  The  effective  modulus  of  first  year  sea  ice  can  be 
related  to  the  brine  content,  Vb  [Timco  and  Weeks,  2010]: 

E  =  10  -  0.3Slvb 

Alternative  formulas  can  be  found  in  Squire  and  Montiel,  [2015]  and  Appendix  A.  If  the 
wave  heights  are  high  enough,  the  bending  stress  will  exceed  the  fractural  strength  or  the 
strain  will  exceed  some  critical  strain.  Unfortunately  there  are  no  definitive 
measurements  of  the  mechanical  properties  of  ice,  and  these  mechanical  properties  are 
function  of  ice  thickness,  age,  porosity,  and  salinity.  Salinity  may  vary  in  the  vertical  cross 
section,  giving  a  single  slab  mechanical  properties  with  depend  on  the  vertical  dimension 
[Squire  and  Montiel,  2015],  For  completeness,  there  is  evidence  that  of  ice  fatigue  due  to 
repeat  loading  may  be  an  important  process  [Langhorne  et  al.,  1998], 

2.2.6  Advection 

On  larger  scales,  say  time  and  space  much  larger  than  typical  wave  periods  and  lengths, 
the  entire  ice  field  may  migrate  as  a  result  of  a  passing  storm.  Brostrom  and  Christensen 
[2008]  (citing  Perrie  and  Hu  [1996])  report  a  migration  up  to  50  km  per  day.  This  amounts 
to  surface  velocity  of  0(0.5  ms'1).  This  migration  is  driven  by  winds  and  currents,  but  wave 
induced  Stokes  drift  and  radiation  stress  must  also  play  a  direct  role.  On  similar  time 
scales,  these  same  forces  may  cause  the  ice  to  converge  and  compress  or  diverge 
depending  on  if  the  forcing  is  directed  into  or  out  from  the  center  of  the  pack  ice, 
respectively.  Of  course,  large  wave  effects  are  only  possible  with  a  sufficient  fetch  which 
necessarily  implies  the  onto  ice  case.  Small,  fetch  limited  waves  are  theorized  to  play  a 
role  in  ice  band  formation  in  the  off-ice  case  [Wadhams,  2000], 

2.3  FEEDBACK 

Waves  may  change  the  nature  of  ice,  and  ice  may  then  change  the  behavior  of  waves.  This 
can  occur  on  very  fine  scales,  say  fracturing  a  single  ice  floe  which  in  turn  modifies  the 
scattering  and  attenuation  the  next  wave  encountered.  On  a  larger  scale,  say  during  the 
passage  of  a  single  storm,  waves  propagate  into  the  MIZ.  Eventually  the  waves  are 
attenuated  by  the  ice,  but  the  waves  are  also  breaking  up  the  ice  which  diminish  one 
mechanism  for  attenuation.  The  waves  continue  to  penetrate  deep  into  the  MIZ  by 
fatiguing  the  ice  along  the  way  until  (1)  the  ice  is  so  strong  that  the  waves  may  no  longer 
pass  or  (2)  the  storm  ends.  This  is  a  sort  of  equilibrium  balance  based  on  the  strength  of 
the  ice  and  the  power  of  the  waves  which  enter  into  the  MIZ.  In  this  way  the  MIZ  is  always 
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in  a  state  of  non-equilibrium,  changing  based  on  the  local  environmental  conditions 
[Wadhams,  2000], 

This  situation  sets  up  the  possibility  of  feedback  loop  in  the  system  on  a  storm  to  storm 
scale.  Continuing  with  the  example  of  a  storm  passing  through  the  MIZ  from  the  open 
water.  Waves  generated  during  the  storm  will  alter  the  MIZ.  If  through  wave  forcing  there 
is  significant  melting/transport  of  ice  would  increases  the  fetch.  An  expanded  fetch 
increases  the  potential  of  more  energetic  wave  event  occurring.  If  another,  very  similar 
storm  follows  in  the  same  area,  the  waves  will  be  bigger  and  the  change  in  the  ice  that 
more  dramatic,  and  the  possibility  of  a  feedback  loop  seems  very  serious  [Asplin  et  al., 
2012;  Asplin  et  al.,  2014;  Thomson  and  Rogers,  2014;  Collins  et  al.,  2015],  Work  is  on¬ 
going  to  determine  the  importance  of  this  feedback  loop  in  the  Arctic  Ocean. 


34 


3  WAVE  EVENT 

3.1  SVALBARD  2010 

The  R/V  Lance  was  in  the  vicinity  of  Hopen  Island  (Svalbard,  NO)  deep  into  the  pack  ice 
when  it  experienced  a  remarkable  wave  event.  A  full  description  of  the  wave  event  and 
measurements  may  be  found  in  Collins  et  al.  [2015],  The  major  findings  of  the  study  was 
an  observation  of  binary  behavior  of  wave  -  ice  interaction  for  the  peak  waves.  Meaning 
that  the  observations  indicated  that  the  peak  waves  were  either  not  allowed  (zero  peak 
wave  energy  transmission)  while  the  ice  was  unbroken,  then  quickly  transitioned  to  total 
transmission  of  peak  wave  energy  with  fractured  ice  floes.  Although  low  frequency  waves 
were  allowed,  the  high  frequency  (but  low  energy)  waves  were  filtered  according  to  the 
local  floe  size  distribution.  The  change  in  local  floe  size  distribution  was  a  result  of 
inhomogeneous  (spatial  gradient)  and  non-stationary  (wave  induced  ice  fracturing  over 
time)  ice  conditions,  though  the  relative  contributions  could  not  be  determined. 

Here  we  take  up  a  few  aspects  of  the  wave  event  which  were  not  fully  explored  in  Collins 
et  al.  [2015],  The  first  is  an  estimation  of  detailed  ice  properties  before  the  fracture  event. 
This  allows  for  an  estimation  of  the  flexural  strength  of  the  ice  to  make  some  inference 
about  the  necessary  wave  conditions.  The  second  is  the  estimation  of  attenuation  after 
the  break  up  event. 

3.1.1  Speculative  Fracture  Investigation  (A.  Marchenko) 

Elastic  ice  is  broken  by  waves  when  tensile  stress  at  the  ice  surface  or  bottom  reaches 
critical  value,  <jcr,  which  is  2-3  times  higher  the  flexural  strength,  <jf  [Marchenko  et  al., 
2014],  Maximal  stress  at  the  surface  of  elastic  ice  sheet  is  calculated  with  the  formula  cr  = 
O.SEhAk2,  where  E  is  the  effective  elastic  modulus,  h  is  the  ice  thickness,  A  is  the  wave 
amplitude,  and  k  is  the  wave  number.  The  effective  elastic  modulus  is  estimated  with  the 
formula  of  Vaudrey  [1977],  E  =  5.31  GPa  —  0A36^v^  GPa,  where  Vb  is  the  liquid  brine 
content  of  the  ice  calculated  with  Frankenstein  and  Garner  [1967]  formula  Vb  = 
S(48.185/|T|  +  0.532).  Mean  values  of  sea  ice  temperature,  T,  salinity,  S,  and  flexural 
strength  measured  in  the  Van  Mijen  Fjord  (Spitsbergen)  in  1-3  May  2010  were  T  =  -1.93°C, 
S  =  3.9  ppt,  and  cy=142.5  kPa.  Therefore  liquid  brine  content  was  Vb  =  100  ppt,  and  E  = 
0.95  GPa.  Critical  tensile  stress  at  the  ice  surface  is  estimated  as  285  <  <jcr  <  427.5  kPa.  The 
observed  wave  length  was  125.6  m  corresponding  to  the  wave  number  k  =  0.05  nr1. 
Assuming  ice  thickness  h  =  0.6  m  the  critical  wave  amplitude  is  estimated  as  0.37  <  A  < 
0.56  m.  Waves  with  higher  amplitudes  should  fracture  sea  ice  immediately.  Langhorne  et 
al.  [1998]  investigated  the  reduction  of  flexural  strength  on  60%  due  to  the  fatigue  failure 
effect.  In  our  case  it  means  a  reduction  of  critical  wave  amplitude  to  the  range  0.22  <  A  < 
0.34  m.  Waves  with  higher  amplitudes  will  fracture  sea  ice  by  low  frequency  fatigue 
during  10  to  100  cycles  of  the  wave  action.  Wave  period  in  our  case  is  12.5  s.  Therefore 
representative  time  of  the  fatigue  fracture  is  varying  from  125  s  to  20  min. 
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3.1.2  Post  Fracture  Attenuation  Rates 

Attenuation  rates  are  calculated  as  a  function  of  period,  T,  using  the  following 
relationship: 

^  ln(Sn+1/Sn) 

«nCO  = - t: - 

Here  Sn+i{/)  is  the  attenuated  wave  spectra  of  original  form  Sn{/).  S(f)  is  either  measured  of 
provided  by  model  application.  Dn  is  the  effective  distance  through  ice  which  the  spectra 
has  propagated.  A  major  assumption  is  stationarity  of  the  wave  and  ice  fields  during  the 
analysis  period  (here  1  hour).  Below  we  explore  different  methods  of  determining  Dn  and 
use  both  measured  and  modeled  S(f). 


3. 1.2.1  METHOD  1:  ALONG  THE  SHIP  TRACK 

Starting  with  the  four  spectra  measured  below.  Si,  S2,  S3,  and  S4  correspond  to  spectra 
measured  from  R/V  Lance  at  02  23:24,  03  00:30,  03  02:30,  and  03  03:30,  respectively. 


Spectral  density 
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Figure  8  Spectral  evolution  over  the  course  of  five  hours.  The  black  line  is  the  SWAN  "no  ice"  reference. 
The  colored  spectra  are  calculated  from  the  sea  surface  elevation  signal  from  R/V  Lance  and  move 
forward  in  time  from  dark  blue  to  green. 

Where  in  the  attenuation  equation,  Dn  is  the  distance  between  the  mean  position  during 
measurement  of  Sn  and  the  mean  position  during  measurement  of  Sn+i.  This  assumes  that 
the  waves  are  traveling  in  the  opposite  direction  of  the  vessel.  There  is  not  a  convincing 
argument  justify  this  assumption,  however  it  should  result  in  at  least  the  right  order  of 
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magnitude.  For  our  three  cases  D\  -  10.1  km,  D2  =  18.7  km,  D3  =  9.9  km.  Leading  to  the 
following  attenuation  rates  as  a  function  of  wave  period  (T). 


Period  (T)  [s] 

Figure  9  The  spectral  attenuation  rates  as  a  function  of  wave  period.  The  light  colored  lines  are  the 
individual  calculations  and  the  black  line  is  their  mean.  The  dotted  line  is  the  empirical  model  from 
Meylan  et  al.  [2014]. 

Shown  is  the  three  attenuation  rates  calculated  from  the  spectra  above  in  light  colors  and 
their  mean  shown  with  the  solid  black  line.  The  dotted  black  line  is  Eq.  3  from  Meylan  et 
al.  [2014]: 

,  ,  a  b 
aM{T)=  —  +  — 

Where  coefficients  a  (2.12  x  10"3)  and  b  (4.59  x  10"2)  were  fit  from  their  measurements. 
The  agreement  between  the  mean  measured  a  and  ctm  is  suspiciously  good. 

3. 1.2.2  METHOD  2:  ICE  CONTOURS  FROM  PIPS 

Let  us  consider  an  alternative  definition  of  Dn  from  the  relative  distance  traveled  into  the 
ice.  This  is  done  by  comparing  the  length  of  rays  parallel  to  the  mean  direction  at  the  peak 
of  the  spectrum  as  output  by  WWIII.  This  assumes  that  the  wave  direction  is  the  same  at 
all  frequencies  and  does  not  change  with  entering  the  ice.  We  track  the  change  in  distance 
to  the  75  %  ice  concentration  contour  given  by  an  ice  model  called  PIPS  as  shown  in  the 
figure  below.  The  assumption  here  is  that  the  75  %  ice  contour,  given  in  the  figure  below, 
from  before  the  storm  is  characteristic  of  the  ice  concentration  during  the  storm 
maintains  shape  over  the  course  of  these  4-5  hours. 
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Figure  10  Map  showing  southern  tip  of  Svalbard,  Norway  including  Hopen  Island  (all  marked  in  green). 
The  ice  concentration  is  a  product  of  the  PIPS  model  and  given  by  the  colored  contours  with  the  white 
indicating  the  75  %  concentration  contour.  The  ship  track  is  marked  in  color  and  the  distances  to  the  75  % 
contour  are  given  parallel  to  the  dominant  wave  direction  as  given  by  WW3  and  represented  by  zi,  zi,  zi, 
and  Z4. 

The  length  of  black  line  segment  which  spans  Hopen  Island  is  approximately  34  km  giving 
a  conversion  rate  from  image  pixels  to  km.  This  segment  is  rotated  to  be  parallel  to  the 
white  arrow  showing  wave  direction  while  retaining  the  length.  This  is  then  copied  until  it 
crosses  the  75  %  contour  and  the  remainder  is  measured  in  pixels  and  converted  to  km. 
The  absolute  distances  for  the  four  segments  are  zi=  89.9  km,  Zi=  77.9  km,  zi=  59.4  km, 
and  Z4=  29.3  km.  The  relative  distances  for  our  three  cases  becomes  D\  =  12.0  km,  D2  = 
18.4  km,  D3  =  10.3  km  (compare  to  path  defined  D\  -  10.1  km,  D2  =  18.7  km,  D3  =  9.9  km). 
This  doesn't  significantly  alter  the  attenuation  rate,  but  the  mean  attenuation  rate  fits  the 
empirical  model  slightly  better. 

3. 1.2.3  METHOD  3:  ICE  CONTOURS  AND  SWAN 

Considering  yet  another  way  to  estimate  attenuation  rates  is  to  use  the  spectra  from 
SWAN  (black  spectra  in  Figure  8).  This  will  give  us  four  attenuation  rates  (1  per  measured 
spectra)  vs.  the  three  from  measured  spectra  only.  This  assumption  is  the  SWAN  provides 
an  estimate  of  spectra  from  the  ice  free  open  waters.  The  problem  is  again  to  determine 
an  appropriate  Dn.,  but  the  nuance  here  is  that  relative  distance  is  no  longer  sufficient,  we 
need  to  estimate  absolute  distance  the  waves  have  traveled  into  the  ice.  We  start  by 
using  the  values  of  zn  as  defined  above  as  first  guess. 
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As  simple  way  to  vary  the  distance  Dn  is  to  scale  the  attenuation  rate.  So  that  multiplying 
the  attenuation  rate  by  2  is  effectively  cutting  the  distance,  Dn,  in  half.  Likewise 
multiplying  by  4  is  quartering  the  distance. 


Period  (T)  [s] 

Figure  11  same  as  Figure  9  with  the  addition  of  colored  dashed  lines  representing  the  mean  attenuation 
rates  using  SWAN.  The  colors  blue,  green,  and  cyan  represent  scaling  of  1,  2,  and  4,  respectively. 

Figure  11  shows  the  mean  of  the  4  attenuation  rate  estimates  using  SWAN  spectra.  The 
initial  guess  at  the  Dn  give  attenuation  rates  which  fell  under  the  empirical  model  from  5  - 
15  s  (expected  effective  range)  but  matched  from  15  -  20  s  (outside  the  expected 
effective  range).  By  quartering  the  distance,  we  get  attenuation  rates  that  better  match 
the  empirical  model  from  4  -  8  s  but  slightly  underestimate  the  model  from  8  -  12  s  and 
slightly  overestimate  thereafter.  This  may  indicate:  1)  a  deficiency  in  the  attenuation 
model,  2)  Dn  was  actually  much  shorter  than  indicated  by  PIPS  (possible  mid-way  through 
the  storm),  3)  deficiency  in  SWAN  or  measured  spectra,  and  4)  all  of  the  above. 

For  sake  of  completeness,  the  ice  edge  could  have  been  defined  at  any  contour  (i.e.  50  %, 
20  %)  or  by  scaling  it,  e.g.  integrating  for  distance  X  but  instead  of  using  the  entire 
distance  use  dx*C  where  C  is  concentration  as  a  fraction.  All  of  these  would  have 
increased  the  distance  Dn  and  lessened  the  agreement  in  the  range  of  periods  where 
attenuation  is  expected. 

Calculating  distance  from  the  ice  edge  is  one  part  of  the  problem,  but  the  other  part  is  ice 
is  changing  in  time,  both  due  to  bulk  migration  of  ice  and  the  advection  of  the  ship 
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through  the  ice,  so  it's  possible  that  we  would  get  a  similar  migration  (in  frequency  space) 
of  the  low-pass  filter  even  with  the  ship  staying  in  one  place. 

4  SPECULATIONS  AND  CONCLUSIONS 

4.1  LARGE  WAVE  EVENTS 

Recently,  a  hypothesis  has  been  put  forward  that  large  waves,  in  particular  large  swell, 
behave  quite  differently  in  ice  than  lower  height  swell  and  even  large  wind  sea.  This  is  still 
very  much  new  territory,  but  we  may  make  a  few  speculations  here  in  the  hopes  that 
future  studies  may  be  able  to  confirm  or  deny. 

Several  studies  report  ice  breakup  due  to  wave  forcing  deep  into  the  pack  ice.  The 
account  of  chief  scientist  of  the  R/V  Polarstern,  E.  Augstein,  was  reported  in  Liu  and 
Mollo-Christensen  [1988]  and  further  expanded  on  in  Brostrom  and  Christensen  [2008], 
Approximately  560  km  from  the  ice  edge,  waves  of  1  m  height  and  18  s  period  broke  up 
ice  and  caused  significant  rafting.  The  wavelength  before  the  ice  broke  was  reported  to  be 
smaller  than  the  wavelength  in  broken  ice  which  motivated  his  formulation  of  the 
dispersion  relation  with  compressive  stress. 

Kohout  et  al.  [2014]  claim  that,  for  particular  wave  periods,  wave  heights  greater  than  3  m 
are  linearly  attenuated  by  ice  vs  the  typical  exponential  attenuation  which  they  found  still 
described  waves  of  heights  less  than  3  m.  In  this  and  more  recent  studies  [Kohout  et  al., 
2015;  Li  et  al.,  2015]  they  argued  that  the  difference  in  attenuation  is  due  to  non-linear 
wave-wave  interaction.  That  effectively,  the  lower  frequencies  continue  to  receive  energy 
from  higher  frequencies  as  a  result  of  the  energy  flux  dictated  by  Sni. 

In  our  estimation,  Sni  acts  over  long  time  and  space  scales,  where,  presumably  the  nature 
of  the  ice  will  be  changing  significantly,  so  in  practice  it  is  very  difficult  to  isolate  the 
influence  of  Sni.  Furthermore,  Sni  is  very  weak  for  typical  ocean  swell  (e.g.  Young  and  Van 
Vledder  [1993])  meaning  for  the  classic  case  of  swell  interning  ice,  Sni  is  likely 
unimportant.  There  is,  however,  an  open  question  to  how  a  change  in  dispersion  might 
affect  Sni  (e.g.  Polnikov  and  Lavrenov  [2007]). 

In  contrast  to  the  longer  scale  non-linearity  associated  with  resonant  wave-wave 
interaction,  Sni,  Collins  et  al.  [2015]  brings  attention  to  shorter  scale,  off-resonant  wave- 
wave  interaction  (i.e.  modulation  instability)  first  discussed  in  the  context  of  ice  by  Liu  et 
al.  [1988],  In  terms  of  attenuation,  after  pack-ice  was  broken  into  floes,  there  was  little  if 
any  attenuation  for  large,  low  frequency  waves  (which  is  also  in  agreement  with  results  of 
Doble  and  Bidlot  [2013]).  They  used  Ml  to  explain  how  waves  are  able  to  march  into  pack 
ice  to  begin  with,  but  not  to  explain  the  different  regimes  of  wave-ice  interaction  (i.e. 
differences  in  wave  attenuation).  Low  frequency  waves  of  the  size  measured  in  Collins  et 
al.  [2015]  have  never  before  been  measured  in  the  Arctic  (see  Appendix  B).  The  regimes  of 
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wave  ice  interaction  were  seemingly  controlled  by  the  ice,  i.e.  whether  the  ice  was  solid  or 
fractured.  It  is  possible  that  the  analysis  of  Kohout  et  al.  [2014]  rolls  the  two  together  and 
gets  significantly  different  attenuation  for  large  waves  because  the  ice  is  fractured 
(fracturing  perhaps  requiring  such  large  waves)  effectively  convoluting  the  two  issues  (ice 
state  and  sea  state).  Unfortunately,  there  is  not  enough  experimental  information 
necessary  to  determine  if  there  are  fundamentally  different  interaction  regimes  for  large 
waves  then  for  small  waves,  but  certainly  the  above  studies  suggest  a  real  need  for  more 
comprehensive  observations.  There  are  several  issues  which  need  attention:  (1)  accurate 
measurement  of  the  attenuation/dispersion  during  and  after  the  ice  breakup,  (2) 
modification  of  regular  wave  growth  and  decay  in  ice,  and  (3)  better  observations  of  the 
ice  break  up  process. 

4.2  WAVES  AS  A  SOURCE  OF  HEAT? 

It  has  long  been  known  that  energy  lost  to  waves  from  breaking  goes  primarily  to 
turbulence  through  the  energy  cascade  and  ultimately  to  heat  via  viscosity.  This  heat 
generation  was  recently  estimated  in  the  context  of  surf  zones  by  Sinnett  and  Feddersen 
[2014],  Two  factors  make  it  plausible  that  this  mechanism  contributes  significantly  to  the 
water  temperature  in  the  surf  zone  at  some  beaches.  First,  in  the  surf  zone,  dissipation 
occurs  over  small  horizontal  distances,  so  dissipation  per  unit  area  is  much  larger  than  in 
the  open  ocean,  where  active  whitecaps  cover  only  a  very  small  fraction  of  the  surface  at 
any  instant.  Second,  the  water  is  shallow,  so  the  heat  cannot  be  mixed  into  lower,  cooler 
water.  Neither  of  these  two  factors  necessarily  applies  to  the  MIZ,  where  dissipation 
occurs  over  larger  horizontal  distances,  and  water  is  deeper.  Flowever,  it  offers  an 
intriguing  possibility,  that  there  is  this  additional  mechanism  by  which  waves  contribute  to 
accelerated  ice  melt  or  decelerated  ice  freezing.  In  any  case,  it  should  not  be  totally 
dismissed  and  it  may  be  worthwhile  in  the  future  to  estimate  this  effect. 

Perhaps  more  interesting  is  the  possibility  that  the  repeated  bending  of  ice  sheets  and 
large  ice  floes  by  wave  action  causes  fatigue  (i.e.  hysteresis  or  internal  friction)  in  the  ice, 
which  results  in  heat.  What  effect  would  this  have  on  the  strength  of  the  ice  and  its 
capacity  to  support  further  strain?  Flowever,  one  could  easily  make  the  counter-argument 
that  the  thermal  and  mechanical  effects  of  the  bending  on  the  ice  would  be  so  difficult  to 
estimate  separately  in  practice,  that  they  would  logically  be  lumped  together  in  any 
predictive  model. 

4.3  FINAL  THOUGHTS 

In  this  document  we  have  summarized  the  nature  of  ice  effects  on  waves  and  waves  on 
ice.  From  a  theoretical  point  of  view,  the  ice  effects  on  waves  is  much  more  developed. 

We  explored  several  formulations  of  the  in  ice  dispersion  relation  in  some  detail.  We  used 
an  example  visco-elastic  model  to  better  understand  the  change  in  dispersion  and  the 
expected  attenuation  once  waves  inter  icy  waters.  The  wave  effects  on  ice  are  intuitive: 
the  may  break  ice  up  and  move  it  around,  but  models  of  these  process  are  very  primitive. 
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On  the  measurements  side,  attenuation  is  well  documented  but  there  are  large 
differences  between  studies  due  in  part  to  the  variable  nature  of  ice.  The  nature  of  ice,  its 
physical  and  mechanical  properties,  are  a  stumbling  block  that  needs  to  be  overcome  in 
future  studies.  Measurements  of  dispersion  are  lacking,  but  the  potential  for  making  such 
measurements  with  marine  radar  is  seemingly  good.  Primitive,  two-way  wave  and  ice  floe 
size  distribution  models  are  immerging  which  are  hopefully  the  vanguards  of  future 
operational  systems.  Before  this  can  happen,  two  broad  issues  need  to  be  addressed:  1)  a 
consistent  definition  of  FSD  needs  to  be  determined  and  used  in  the  literature  and  2) 
continued  observations  of  FSD  which  might  be  used  to  compare  to  the  models.  In 
particular  a  description  of  FSD  is  needed  which  either  definitively  supports  or  challenges 
and  supersedes  the  current  power  law  based  models. 

The  wave-ice  interaction  problem  is  difficult,  in  some  cases  intractable,  because  of  the 
disparate  scales  involved  and  the  varied  nature  of  ice  in  the  field.  With  the  ONR  Sea  State 
field  experiment  underway,  there  is  the  looming  prospect  of  high  quality  field  data  and 
new  theoretical  tools  which  we  hope  will  usher  in  new  breakthroughs  in  wave-ice  physics. 
Flere  we  must  face  that  no  matter  how  comprehensive  a  measurement  campaign, 
developing  a  model  which  can  predict  two-way  coupled  behavior  of  waves  and  ice  over  all 
the  various  scales  is  a  daunting  task.  Nonetheless,  we  have  never  been  so  well  equipped 
to  undertake  such  challenges  and  the  prospects  for  progress  are  favorable. 
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APPEDDIX  A:  BENDING  AND  DISPERSION  OF  VISCOELASTIC 
WAVES  (MARCHENKO) 

A.l  BENDING  MOMENT  IN  ICE  PLATE  WITH  ELASTIC  AND  CREEP 
RHEOLOGICAL  PROPERTIES 

Two  dimensional  bending  deformations  of  floating  solid  ice  caused  by  wave  induced  water 
pressure  below  the  ice  are  described  by  equations  of  static  equilibrium 


da  da  .  da  da 

— —  +  .  **  =  0 ,  ^  *z  +  =  p:g , 


dx  dz 


dx  dz 


(la) 


where  x  and  z  are  the  horizontal  and  vertical  coordinates,  cr  ( i,j  =  x,z )  are  the  ice 
stresses,  and  pt  is  the  ice  density.  Ice  inertia  is  ignored  in  comparison  with  inertia  of  the 
water  involved  in  the  wave  motion.  It  is  valid  when  wave  lengths  and  water  depth  are  much 
greater  the  ice  thickness. 

Linearized  boundary  conditions  have  the  following  form 

=~Pa>  =  0 '  z  =  h/2;  a,z  =  pw ,  ax_  =0,  z  =  -hl 2,  (2a) 


where  pa  is  the  atmosphere  pressure,  and  pw  is  the  water  pressure  below  the  ice. 


Figure  A-l.  Ice  plate  floating  on  the  water  surface. 

Multiplying  the  first  equation  (l)on  z  and  integrating  it  over  z  from  z  =  -ht 2  to  z  =  ht 2 
we  find 


8MXX 

dx 
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Integration  of  the  second  equation  (1)  over  z  lead  to  the  equation 

dN x  ,  .  . 

=  Pigh  -  Ap ,  Ap  =  pw  -  pa . 
dx 

Combination  of  equations  (3)  and  (4)  gives 


(3a) 


(4a) 
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=  Pa~  PW 


(5a) 


9Xx 

dx 2 


Link  between  bending  moment  Mxx  and  water  surface  elevation  77  is  established  by  the 
use  of  rheological  equations.  We  assume  plane  strain  state  in  the  ice  with  respect  to  the 
axis  y  perpendicular  to  the  plane  (x,z).  We  assume  creep  rheology  of  ice  for  the  stress 
deviator  and  elastic  rheology  for  the  pressure 

Q( £  — £  ^ 

4  +  4  = - ”  °  -  =  24  +  ■  (M 

at 


where  X  and  //  are  effective  elastic  constant,  and  sel  and  scr  are  elastic  and  creep  strain 
rates.  Elastic  strain  rate  is  determined  by  the  formula 


1  d{(jxx-cr:z) 
2/7  dt 


(7a) 


Elastic  constant  are  expressed  by  the  formulas 


E  Ev 

2(1  +  v)  '  2(l  +  v)(l-2v)  ’ 


(8a) 


where  E  is  the  effective  elastic  modulus,  and  v«0.33  is  the  Poisson's  ratio  of  ice.  The 
elastic  modulus  is  calculated  with  the  formula  (Vaundray,  1977) 

E  =  5.31-0.436^7-  (9a) 


where  vh  is  calculated  in  ppt.  Liquid  brine  content  vb  is  calculated  as  follows  (Frankenstein 
and  Garner,  1967) 


vh  =  Ssi 


48.185 
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,  -22.9°C<r<-0.5°C, 


(10a) 
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where  77  and  Ssi  are  the  temperature  and  salinity  of  sea  ice  averaged  over  ice  thickness. 
Creep  strain  rate  is  expressed  by  the  formula  (Durham  et  al,  1997) 


K=K  I  axx  -  ^ZZ  r  (O'xx  -  0"zz)' 


(11a) 


where  creep  constant  K  could  expressed  in  terms  of  representative  strain  rate  (s0)  and 
stress  (cr0)  in  creep  tests:  K  =  £qct^  .  Assuming  that  ^0=10_6s_1  and  cr0  =  1  MPa  (Cole  et 
al,  1998)  we  estimate  K  =  10“6s-1MPa"4 . 
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It  is  assumed  that  bending  deformations  of  floating  ice  are  caused  by  surface  waves  with 
period  about  10  s  and  wavelength  125  m.  Maximal  elastic  stress  caused  by  wave  induced 
bending  deformations  of  floating  elastic  ice  plate  is  estimated  by  the  formula 


cr 


el 


Ehak2 

2(l-v2)' 


(12a) 


where  a  and  k  are  wave  amplitude  and  wave  number.  Assuming  h  =  0.5  m,  E  =  2GPa  , 
a  =  0.5  m  and  k  =  0.05  irf1  we  find  <jel  =1.25  MPa . 


For  the  comparison  of  the  terms  in  the  left  side  of  the  first  equation  (6)  we  assume  that 
axx  -  a„  =  ael .  The  ratio  of  elastic  strain  rate  to  creep  strain  rate  is  equal 


4  = 

4  EKct],' 


(13a) 


where  &>  =  0.6s-1  is  the  wave  frequency.  Numerical  estimates  show  that  selsj  ~  1.5 -102, 
i.e.  creep  strain  rates  are  much  smaller  than  elastic  strain  rates. 

From  equations  (3)  and  (5)  follows  that  M  ~  a  h2  and  M  k2  ~  p  when  the  ice  is 
deformed  by  a  wave  with  wavelength  k.  Therefore  pw~axxh2k 2  and  pH,«crxx  when 

h2k2  «  1 .  Further  this  fact  is  used  for  the  construction  of  solution  of  rheological  equations 

(6). 

Solution  of  equations  (6)  is  constructed  by  the  method  of  successive  approximations  with 
the  first  term  describing  pure  elastic  solution 


(14a) 
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&o,xx  +  CTo.zz  =  2(A  +  p)(sxx  +  sz:)  . 


Assuming  that  <yQxx  »cr0__  we  find  <r0  YT  =  E(l-v2)  lexx .  Taking  into  account  the 

Kirchoff's  hypothesis  (Timoshenko  and  Woinovsky-Krieger,  1959)  the  final  form  of  the  first 
order  approximation  is 
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where  tj  is  the  deflection  of  the  mid-surface  from  the  horizontal  plane  (Fig.  1). 


In  the  next  order  of  the  approximation  we  assume  that  crxx  =  cr0xx  +crlxv,  where  stress  crlxx 
accounts  the  influence  of  creep  effect  and  calculated  from  the  equation 
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From  equations  (15)  and  (16)  follows  that  bending  moment  is  determined  by  the  equation 
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A. 2  DISPERSION  EQUATION  FOR  THE  WAVES  PROPAGATING  IN  IDEAL  FLUID 
BENEATH  AN  ICE  PLATE  WITH  ELASTIC  AND  CREEP  PROPERTIES 
Linearized  kinematic  and  dynamic  boundary  conditions  at  the  surface  of  ideal  fluid  covered 
by  the  ice  plate  are  written  as  follows 


dr]  _  dcp  dcp 
dt  dz  '  dt 


+  gf 1 


1  d2Mxx 
p  dx2 


=  0,  z-0, 


(18a) 


where  cp  is  the  velocity  potential  in  the  fluid,  i]  is  the  water  surface  elevation,  and  p  is  the 
ice  density.  Bending  moment  Mxx  is  expressed  by  i]  using  equation  (17).  Boundary 
conditions  (18)  are  reduced  to  the  boundary  conditions  at  the  surface  of  ideal  fluid  covered 
by  elastic  plate  when  Dcr  =  0  . 


Substituting  in  equations  (17)  standard  expressions  for  the  water  velocity  potential  and 
water  surface  elevation  describing  surface  waves  in  the  fluid  with  depth  H 


^^cosl m^_M  +  c.C., 


cosh[AW] 


r]  =  Ad 


i{kx+cot) 


+  C.C., 


(19a) 


we  find  the  dispersion  equation 

c o 2  =  (o2fgw  +  i\A |4A ,  A  =  Dcrks/p,  (20a) 

(o2fgw  =  k  tanh[A-//](g  +  Delk4  Ip). 
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APPENDEX  B:  WAVE  HEIGHTS  MEASURED  IN  ARTIC  ICE 

Here,  to  support  our  claim  to  have  measured  the  largest  waves  in  ice  in  the  Arctic  region, 
we  have  attempted  to  make  an  exhaustive  review  of  the  literature.  An  earlier  version  of 
this  appendix  was  published  as  supplemental  material  for  [Collins  et  al.,  2015], 

Studies  report  wave  height  using  a  variety  of  metrics.  In  studies  where  the  individual 
maximum  wave  height,  Hmax,  or  the  root-mean-square  wave  height,  Hrm s,  was  reported, 
we  use  a  simple  formula  to  convert  to  the  spectral  significant  wave  height,  Hm o. 
Considering  typical  wave  periods  and  analysis  times,  the  following  relationships  derived 
from  a  Rayleigh  distribution  are  roughly  applicable  (e.g.  Holthuijsen  [2007]). 

^m0  —  1.4 Hrms  ~  O.SH max 

When  only  the  energy  (or  variance)  density  is  reported  for  a  particular  spectral  band,  an 
attempt  is  made  to  estimate  a  JONSWAP  spectrum  [Hasselmann  et  al.,  1973]  based  on  the 
energy  level  at  the  peak  frequency.  This  is  done  in  a  conservative  manner  such  that  the 
estimation  of  Hmo  is  an  upper  limit  (or  perhaps  slightly  overestimated). 

We  tried  to  be  thorough,  but  invariably  there  will  be  studies  which  were  passed  over  or 
could  not  be  verified.  One  that  fits  into  the  latter  category  is  work  by  P.  Wadhams,  in 
which  the  measurement  of  wave  heights  was  made  by  an  inverted  echo  sounding  from  a 
submarine  [Wadhams,  1972;  Wadhams,  1978],  Though  these  citations  are  indexed  (e.g. 
google  scholar,  web  of  science)  actual  records  (proceedings  in  the  first  case,  an  article  in 
the  second)  could  not  be  located. 

As  far  as  the  Arctic  region  (see  Table  SI)  is  concerned,  our  study  is  by  far  highest  and  is 
similar  to  wave  heights  observed  in  the  ice  free  measurements  of  Thomson  and  Rogers 
[2014]  in  the  Beaufort  Sea.  As  it  may  be  of  interest,  we  have  also  tabulated  the  cases  of 
waves  measured  in  ice  in  the  Antarctic  as  well  (see  Table  S2). 
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Table  2  Summary  of  the  published  accounts  of  waves  measured  in  ice  in  the  Arctic  region. 


PUBLICATION 


STUDY 

YEAR 


LOCATION 


WAVE 

METHOD  HEIGHT 

(Hmo) 


[Fakidov,  1934] 

1933- 

1934 

Arctic  Ocean 

Visual 

«  1  m 

[Crary  et  al.,  1952] 

1951 

Arctic  Ocean 

Gravity  Meter 

5xl0'4  m 

[Hunkins,  1962] 

1957- 

1958 

Beaufort  Sea 

Gravity  Meter6 

«  1  m 

[Wadhams,  1975] 

1972 

Canadian  North 
Atlantic7 

Airborne  Laser 

~2  m 

[Squire  and  Moore, 
1980] 

1979 

Bering  Sea 

floe-borne 

accelerometers 

~1.5  m 

[Wadhams  et  al., 
1986;  Wadhams  et 
al.,  1988] 

1978, 
1979,  & 
1983 

Bering  Sea  and 
Greenland  Sea 

Buoys 

~1  m 

[Liu  et  al.,  1991] 

1989 

Labrador  Sea 

Buoy 

~2.6  m 

[McKenna  and 
Crocker,  1992] 

1989 

Labrador  Sea 

floe-borne 

accelerometers 

~0.5  m 

[Rottier,  1992] 

1989 

Barents  Sea  and 
Fram  Straight 

Buoy 

-0.75  m 

[S.  Frankenstein  et 
al.,  2001]8 

1990 

_ 

Barents  Sea 

floe-borne 

accelerometers 

? 

[Marko,  2003] 

1998 

Sea  of  Okhotsk9 

ADCP 

~1.5  m 

[Asplin  et  al.,  2012] 

2009 

Beaufort  Sea 

3-D  ship-borne 
recorder 

~0.75  m 

This  Study 

2010 

Barents  Sea 

Ship  GPS 

4-5  m 

6There  are  other  studies  which  analyze  the  observations  of  gravity  meters,  seismometers,  or 
tiltmeters  deep  in  the  ice  field  (see  the  references  in  Wadhams  [1975]),  all  of  which  record  small 
(O(lxl0  4  m))  vibrations  of  the  ice. 

7Not  within  the  Arctic  region. 

8An  event  very  similar  to  the  one  reported  in  out  manuscript  was  described,  but  measurements 
were  not  made  during  it. 

9Not  within  the  Arctic  region. 
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Table  3  Summary  of  the  published  accounts  of  waves  measured  in  ice  in  the  Antarctic  region. 


PUBLICATION 


STUDY 

YEAR 


LOCATION 


WAVE 

METHOD  HEIGHT 

( Hmo ) 


[Robin,  1963] 

1959- 

1960 

Wendell  Sea 

Ship-borne  Recorder 

~0.5  m 

[Liu  and  Mollo- 
Christensen,  1988] 

1986 

Wendell  Sea 

Visual 

~1  m 

[Crocker  and 
Wadhams,  1988] 

1986 

Ross  Sea 

Wire  Strainmeter 

«  1  m10 

[Fox  and  Haskell, 
2001] 

1998 

Ross  Sea 

Floe-borne  recorder 

?n 

[Hayes  et  al.,  2007] 

2003 

Bellinghausen  Sea 

ADCPonan  AUV 

~3  m 

[Doble  and  Bidlot, 
2013] 

2000 

Wendell  Sea 

Buoy 

~10  m 

[Meylan  et  al.,  2014; 
Kohout  et  al.,  2014] 

2012 

Antarctic  Ocean 

Floe-borne  recorder 

~6  m 

'“Measured  strain  (not  amplitude),  a  conversion  depended  on  Young's  modulus  (not  measured). 
"Measurements  shown  only  in  terms  of  acceleration  and  acceleration  spectra. 
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